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ABSTRACT 

The  Finite  Element  Analysis  Prograa  (FEAP)  was  modified  and 
integrated  with  the  Naval  Postgraduate  School  version  of  the 
Computer  Analysis  Language  (CAL-NPS) .  This  enables  the 
solution  of  linear  and  non-linear,  two  and  three  diaensional 
heat  conduction  problems  in  an  interactive  node.  The  usual 
types  of  boundary  conditions,  including  radiation,  nay  be 
specified.  The  heat  conduction  group  includes  proapts  for 
user  supplied  data.  Several  existing  CAL-HPS  coaaands  were 
iaproved  and  a  "HELP"  facility  was  added.  Coaaands  were 
added  for  visual  display  of  the  finite  eleaent  aesh  at 
graphics  terainals.  The  User  Guide  for  thA«:  expanded 
version  of  CAL-HPS  is  provided. 
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I.  IHTBODOCriOS 


A.  GSHEHAL  DBSCBIPTIOH 

Since  the  implementation  of  the  Finite  Element  Analysis 
Prograa  (hereafter  refers  d  to  as  FEAP)  at  the  Naval 
Postgraduate  School,  it  has  had  only  limited  use.  it  is  now 
available  to  general  NPS  users  in  an  interactive  fora 
through  the  IBS  3033  7a/HVS  time  sharing  system.  The  data 
management  and  user  interactive  command  structures  were 
established  within  the  existing  interactive  prograa. 
Computer  Analysis  Language  (CAL).  This  integration  of 
systems  provides  the  ability  to  solve  linear  and  nonlinear, 
steady  and  unsteady,  two  and  three  dimensional  heat  conduc¬ 
tion  problems  involving  temperature  dependent  theraophysical 
properties  and  complicated  radiation/convection  boundary 
conditions. 

Additional  capability  was  provided  to  CAL  system  users 
both  in  changes  to  existing  subroutine  groups  and  in  the 
addition  of  a  new  graphics  group  with  its  attendant  command 
structure.  The  graphics  functions  enable  the  user  to  plot 
two-  and  three-  dimensional  structural  and  heat  transfer 
meshes.  Another  new  facility  is  the  HELP  operation  which 
allows  a  user  experiencing  trouble  with  a  particular  opera¬ 
tion  to  interactively  obtain  assistance. 

B.  HISTORICAL  BACKGROUND 

The  original  CAL  program  was  developed  by  Professor  E. 
L.  Hilson  of  the  University  of  California  in  1977  [Hef.  1  ]. 
It  was  later  adapted  and  modified  for  use  at  the  Naval 
Postgraduate  School  by  LCDS  L.  B.  Elliot  2].  The  PEAP 

prograa  was  written  by  Professor  R.  L.  Taylor  of  the 


University  of  California  [fief.  3]  in  1977.  Implementat ion 
of  FEAP  at  Naval  Postgraduate  School  was  done  by  LT  J.  H. 
Bettencourt  [Ref.  *]. 

C.  OBJECTIVES 

The  objectives  of  the  author's  work  have  been  to: 

1.  integrate  the  data  sanagesent  systes  of  FEAP  with  CAL 
to  create  an  interactive  conduction  heat  transfer 
probles  solving  systes; 

2.  modify  existing  operations  to  extend  their  useful¬ 
ness; 

3.  add  operations  to  extend  the  capabilities  of  the 
program  to  include  graphics; 

tt.  create  a  HELP  facility; 

5.  create  a  USER'S  MANUAL  to  facilitate  use  of  this 
program. 
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II.  ORGANIZATION  OF  BEAT  TRANSFER  GBOPP 

This  chapter  provides  a  general  overview  of  the  organi¬ 
zation  of  the  Heat  Transfer  Group  of  CAL.  It  is  intended  to 
provide  sufficient  information  to  permit  users  to  operate 
the  CAL  heat  transfer  package. 

The  execution  of  the  program  is  flexible  and  controlled 
by  user  selection  of  operations  in  a  logical  sequence  from 
the  commands  that  are  available.  There  are  two  broad 
categories  of  operations,  data  input  and  problem  solution. 

A.  DATA  INP0T 

1 .  Initializat  ion 

The  heat  transfer  group  performs  matrix  creation  and 
manipulation  automatically.  As  the  problem  progresses, 
arrays  are  created,  altered  and  deleted  under  program 
control.  Through  the  HTXFR  operation  the  user  provides 
sufficient  information  to  establish  the  initial  arrays  for 
data  input  and  problem  solution.  The  number  of  nodes, 
number  of  elements,  number  of  material  sets,  spatial  dimen¬ 
sion,  number  of  degrees  of  freedom  per  node  and  the  maximum 
number  of  nodes  per  element  are  required.  For  heat  transfer 
problems  the  number  of  degrees  of  freedom  per  node  is  always 
one.  The  option  to  assign  a  higher  number  is  available 
because  the  equation  solvers  in  this  program  are  applicable 
to  other  fields  of  which  future  work  may  make  use. 

2.  Coordinates 

Nodal  coordinates  are  input  via  the  COOBD  operation. 
This  operation  has  built  in  node  generation  capability.  By 
specifying  an  initial  point  and  a  node  generation  vector. 


the  aser  may  easily  input  large  aeshes.  Coordinate  systee 
conversion  is  also  available.  Coordinates  nay  be  input  in 
the  Cartesian  systee,  the  cylindrical  systee  with  any  one  of 
the  three  axes  longitudinal,  the  spherical  systee  or  any 
coebination  of  the  above  systems.  &11  coordinates  are 
converted  to  cartesian  coordinates  for  use  with  CAL. 

3.  Slaisflt  C2aag£&i 111 

The  ELCOH  operation  inputs  the  element  connectivity 
data.  Here  again  is  a  generation  capability.  The  user  nay 
specify  the  connectivity  for  one  eleaent  and  a  generation 
vector  to  create  additional  rows. 

for  the  two-diaensional  elements  the  user  may 
specify  a  4  to  9  node  isoparametric  eleaent.  There  is  an  8 
to  21  node  isoparametric  eleaent  for  three-dimensional 
elements.  Both  of  these  elements  must  follow  the  numbering 
convention  shown  in  Appendix  A. 

9.  flaJLSEiAl  Properties 

The  required  amount  of  material  property  information 
varies  from  problem  to  problem.  The  PROP  operation  prompts 
the  user  for  the  information  required  to  solve  the  problen 
at  hand.  At  a  minimum  the  material's  conductivity  (k) , 
specific  heat  (c)  ,  specific  mass  (rho)  rheat  generation  per 
unit  volume  (q'*')  and  the  geometry  type  (plane  or  axisyae- 
tric)  must  be  specified  in  a  consistent  system  of  units. 
Appendix  A  includes  examples  of  consistent  systems. 

The  user  also  inputs  the  number  of  Saussian  points 
per  direction  for  quadrature  and  codes  for  temperature 
dependent  properties  and  boundary  conditions.  The  codes 
indicate  other  information  required.  Temperature  dependent 
properties  are  input  as  tables  and  linear  interpolation  is 
used  to  determine  the  property  value  at  a  given  temperature. 
Boundary  conditions  are  specified  as  shown  in  Appendix  A. 
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If  a n  exterior  boundary  line  or  surface  condition  is  not 
specified,  it  is  assuaed  to  be  insulated. 

5.  Constant  Teaperature  go^es 

For  probleas  requiring  certain  nodes  to  be  at 
constant  teaperature,  the  CTEMP  operation  is  available. 
This  operation  say  also  be  used  daring  the  solution  stage  of 
the  problea  to  provide  step  changes  at  previously  specified 
constant  teaperature  nodes.  Because  this  operation  gener¬ 
ates  an  array  used  in  profiling  the  solution  equations,  the 
user  aay  not  change  the  node  nuabers  that  were  established 
as  constant  constant  teaperature  nodes  after  execution  of 
the  PROF  operation.  The  teaperatures  of  these  nodes, 
however,  aay  be  changed. 

The  PROF  operation  establishes  the  equation  profile 
for  problea  solution.  Prior  to  the  execution  of  this  opera¬ 
tion,  any  data  input  aay  be  changed  by  specifying  the  appro¬ 
priate  operation  and  re-entering  the  data.  After  its  execu¬ 
tion  the  user  aay  not  change  the  nodes  designated  as 
constant  temperature  nodes  to  teaperature  varying  nodes  or 
vice  versa.  The  user  aay  change  the  value  of  the  constant 
teaperatures. 

B.  SOLOTIOI 

The  aatrix  formulation  of  the  heat  transfer  problea  as 
discussed  in  Reference  4  is: 

[K]{T}  ♦  [C](T}  ♦  {F}  =  CO}  O) 

where  (K)  represents  the  conductivity  aatrix,  (C)  represents 
the  heat  capacity  aatrix,  (F)  represents  the  flux  vector  and 
{T}  represents  the  temperature  vector.  The  derivative  of 
{T}  with  respect  to  time  is  (T)  .  The  flux  vector  includes 


heat  generated  per  unit  volume  and  boundary  fluxes  as  speci¬ 
fied  for  given  boundary  surfaces.  This  is  a  fully  general¬ 
ized  formulation,  including  non-linearities,  since  the 
mat  rices  C  K]  and  [C]  and  the  vector  {P}  can  be  temperature 
depen  dent. 

laoisa  saaia s&iriti  aaiiii 

Por  heat  transfer  problems,  which  involve  only  one 
degree  of  freedoa  per  node,  the  conductivity  matrix,  [K], 
will  always  be  syaaetric.  The  capability  for  generating 
unsyaaetric  matrices  was  provided,  but  will  only  be  appli¬ 
cable  when  additional  types  of  probleas  are  programmed  into 
the  CAL  system.  The  command  for  the  unsymaetric  conduc¬ 
tivity  aatrix  formulation  is  USIMC. 

2 •  Easiing  SS3£  capacity  Matrix 

The  heat  capacity  aatrix,  [C],  used  in  time  depen¬ 
dent  probleas,  can  be  generated  with  either  CCAP  or  LCAP. 
To  fora  a  consistent  capacitance  approximation  use  the  oper¬ 
ation  CCAP.  A  lumped  capacitance  approximation  is  formed 
using  the  LCAP  operation. 

3.  ZSxaiM  Zlai  Vector 

To  complete  the  problem  formulation,  the  flux 
vector,  JP},  must  be  generated.  The  FORM  operation  forms 
the  flux  vector  taking  into  account  the  internal  heat  gener¬ 
ation  and  boundary  surface  fluxes  as  indicated  in  the  PROP 
operation. 

* •  laiidtioa  Soilin'! 

Once  the  time  independent  problem  is  formulated,  the 
temperature  vector,  (T}  ,  is  calculated  by  the  CALC  opera¬ 
tion.  Time  dependent  probleas  do  not  use  this  operation, 
but  rather  the  ordinary  differential  equation  solver. 
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5.  fim  Qtist  Qiliaatz  Differential  gguatjon,  Sglver 

The  first  order  ordinary  differential  equation 
solver  is  accessed  with  the  ODE  operator.  It  eaploys  the 
Zienkiewicz  two-  and  three- level  scheaes  [Ref.  4V  5]. 

In  addition  to  the  tiae  step  size  change  using  the 
DTIH  operation,  an  optional  autoaatic  tiae  step  adjustaent 
is  incorporated  in  the  ODE  operation.  The  nora  of  the 
difference  between  teaperature  vectors  at  two  consecutive 
tiaes  is  coaputed  at  each  step.  If  the  nora  is  less  than  a 
user  specified  aaziaua  teaperature  difference,  the  tine  step 
is  doubled  before  going  to  the  next  step.  If  the  nora  is 
greater  than  a  user  supplied  niniaun  teaperature  difference, 
the  tiae  step  is  halved  and  calculation  for  that  tiae  step 
is  repeated  until  the  nora  is  acceptable.  If  the  tenpera- 
ture  differences  are  specified  as  zero,  no  tiae  step  adjust- 
nent  will  be  perforned. 

The  user  specifies  cne  of  three  functions  which  are 
perforned  by  the  ODE  operation.  A  second  operation  nane, 
which  aust  be  separated  by  a  coaaa ,  follows  ODE.  The 
options  are  INIT,  LINE  or  9 DAD. 

The  operation  ODE, INIT  is  used  to  input  the  integra¬ 
tion  constants  theta,  beta  and  ganma  [Ref.  5],  the  aaxiaua 
and  ainiaua  teaperature  differences  for  the  autoaatic  tiae 
st9 p  adjustaent  and  the  initial  teaperature  vector.  No  tiae 
integration  is  perforaed  by  this  instruction. 

ODE, LINE  performs  the  two  point  schene  and  the 
current  teaperature  vector  is  substituted  by  the  newly 
calculated  teaperature  vector. 

The  ODE, QUAD  operation  is  siail'r  to  the  ODE, LINE 
operation  but  uses  the  three  point  scheae. 

6.  Hinting  Ssdfii  leaperatures 

Once  the  teaperature  vector  has  been  updated,  the 
PTEMP  operation  prints  the  teaperatures  in  node  nuaber 
order. 
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III.  ORGAHIZATIOR  OP  GRAPHICS  GROUP 


The  graphics  group  is  capable  of  displaying  aeshes 
either  on  the  IBS  3277  dual  screen  terainal  system  or  any 
PLOT- 10  compatible  terainal.  It  is  initiated  through  the 
use  of  the  GRAPH  operation  by  which  the  user  specifies  the 
type  of  graphics  terainal  in  use. 

A.  TITLE 

The  TITLE  operation  is  used  to  title  the  aesh  being 
displayed.  The  user  aay  input  up  to  three  lines  of  fifteen 
characters  for  the  title.  This  operation  aust  be  specified 
iaaediately  prior  to  the  displaying  operation. 

This  operation  calls  the  OSRIH  subroutine  which  reads 
froa  the  terainal  three  lines  of  characters.  A  flag  (IFLAG) 
is  set  to  indicate  a  title  is  to  be  printed. 

B.  HEAT  TRANSFER  HESH 

The  PLHX  operation  locates  the  heat  transfer  coordinate 
and  connectivity  arrays  and  displays  the  aesh.  The  viewing 
area  is  optimized  so  the  longest  dimension  is  full  screen. 
The  aaziaua  and  ainiaua  coordinate  values  for  each  direction 
are  displayed.  This  operation  displays  two-  and  three- 
dimensional  aeshes.  If  the  aesh  is  three-dinensional,  the 
user  aust  specify  the  viewing  plane. 

This  operation  calls  either  the  FPPLOT  (2-D)  or  FP3PLT 
(3-D)  subroutine.  These  subroutines  scan  the  coordinate 
arrays  for  aaziaua  and  ainiaua  values  and  initialize  the 
graphics  screens.  They  both  use  subroutine  BOX  to  set  the 
virtual  window,  set  the  screen  window  and  draw  a  box  around 
the  plotting  area.  The  standard  element  connectivity  is 


stored  is  an  array  which  is  used  in  conjunction  with  the 
user's  eleaent  connectivity  matrix  to  draw  the  eleaents  line 
by  line.  The  center  node  is  plotted  using  a  syabol. 
The  user's  title  is  plotted  (if  IFLAG  is  not  0)  by  subrou¬ 
tine  OSRTIT,  which  then  resets  IFLAG.  The  aaxiaua  and 
ainiaua  values  of  the  vertical  and  horizontal  coordinates 
are  printed  by  SC BOAT .  Subroutine  SCB DAT  also  indicates  axis 
orientation.  Prior  to  terainating  the  screen  graphics, 
subroutine  TITLE  is  used  to  write  the  title  box  identifying 
the  type  of  aesh  being  displayed. 

C.  STRUCTURAL  HESH 

The  PLST  operation  displays  the  structural  systea  speci¬ 
fied  by  user  supplied  coordinate  and  connectivity  arrays. 
The  viewing  area  is  optiaized  as  described  for  the  PLHX 
operation.  Likewise,  the  user  nay  specify  the  viewing  plane 
as  either  the  X-Y,  Y-Z  or  X-Z  plane. 

This  operation  calls  the  CLPLOT  subroutine.  CLPLOT  uses 
the  saae  logic  and  subroutines  as  FPPLOT.  The  difference 
between  the  two  subroutines  is  in  the  aethod  of  storage  of 
the  coordinate  arrays  and  the  connectivity  matrix. 
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17.  CHARGES  TO  BUSTING  GROUPS 

Hodifi  cat  ions  were  nie  in  several  existing  subroutine 
groups  to  sake  their  operations  aore  versatile.  One  utility 
subroutine  was  iaproved  and  one  was  added. 

The  changes  that  were  nade  are  sensitive  to  prior 
versions  of  CAL.  The  saae  results  are  obtained  for  previ¬ 
ously  existing  operation  coaaands.  Ho  files  used  with  other 
editions  of  CAL  need  to  be  aodified  to  operate  with  this 
prograa. 

A.  UTILITY  SUBROUTINES 

1.  Subroutine  RCARD 

Subroutine  RCARD  reads  and  interprets  the  operation 
coaaands.  An  operation  coaaand  has  the  fora 
OP, HI , M2, M3 ,M4 ,M5 ,N1, N2,N3,N4 

where  OP  is  the  operation  naae,  HI  to  H5  are  aatrix  naaes 
and  HI  to  H4  are  integers.  Previously,  HI  to  N4  had  to  be 
values  greater  than  or  equal  to  zero.  The  syabol/state 
logic  aatrix  and  subsequent  action  codes  were  aodified  to 
allow  users  to  input  negative  integer  values. 

This  is  iaportant  ,  for  the  ZERO  operation  which 
enables  users  to  create  aa trices  with  a  given  value  in  the 
diagonal  locations  and  another  value  in  the  off  diagonal 
locations.  Prior  to  this  change,  a  user  desiring  to  create 
such  a  aatrix  with  negative  values  had  to  input  the  aatrix 
row  by  row. 

2.  lilfelfllftlag  OISflX 

Subroutine  PRTCHX  is  a  new  utility  subroutine  which 
allows  the  CAL  prograa  to  invoke  most  CP/CMS  coaaands. 
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After  the  invoked  command  is  executed,  control  returns  to 
CAL. 

It  is  presently  used  in  tie  SAVE  and  RESUME  opera¬ 
tions  to  invoke  the  CHS  coeeand 

PIL2DEP  HSAVE  DISK  HI  SAVE  (HECPH  VS  LRECL  7290 
BLRSIZE  7294) 

where  HSAVE  is  a  logical  unit  number  assigned  to  the  SAVE 
and  RESUME  operations  and  Si  is  a  user  input  file  nane. 

B.  SBOOP  1 

1 •  §1 IE  Q£S£fttjga 

The  save  operation  creates  a  file  on  the  user*s 
A-disk  containing  all  arrays  in  storage  at  the  tiae  of 
issuance.  Previously  the  entire  100,000  word  aain  array  was 
stored,  regardless  of  how  aany  words  were  actually  being 
used.  It  was  always  stored  under  the  naae  FILE  02, 

preventing  the  user  from  saving  more  than  one  problem. 

The  net  hod  of  storing  the  array  was  changed  to  store 
just  the  number  of  locations  actually  being  used.  More  than 
one  problem  may  be  saved  because  the  subroutine  FRTCMX  was 
used  to  create  a  SAVE  file  with  a  naae  assigned  by  the  user. 
If  a  name  is  not  specified  PILE  02  will  be  the  name  of  the 
saved  problem. 

2 •  ££§2fl£  operat ion 

The  RESUME  operation  reads  a  saved  file  into  memory. 
It  was  altered  to  read  named  files  saved  by  the  new  SAVE 
operation.  If  a  naae  is  not  specified  FILE  02  will  be  read. 

3.  LOADI  operation 

The  LOADI  operation  loads  integer  arrays.  The  arays 
were  input  row  by  row.  The  option  to  generate  arrays  was 
added.  The  user  may  specify  one  row  and  a  row  generation 
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vector.  The  number  of  rows  specified  will  be  automatically 
generated.  This  operation  was  aoved  froa  the  static  anal* 
ysis  subroutine  group  to  the  general  aatriz  coaaand  group. 

«.  s&isi  aagcatiaa 

The  PRINT  operation  prints  an  array  in  aatriz 
foraat.  Previously  it  could  only  print  arrays  containing 
real  numbers,  resulting  in  the  erroneous  printing  of  arrays 
containing  integers.  The  user  aay  now  specifiy  whether  the 
array  to  be  printed  contains  real  or  integer  values. 

5-  fljil  Operation 

The  HELP  operation  was  added  to  the  general  coaaand  group 
and  provides  the  user  with  information  on  the  use  of  all  the 
available  operations.  It  accesses  a  file  of  instructions, 
sorts  through  thea  to  find  the  desired  operation  and 
displays  the  appropriate  information  on  the  screen. 

C.  GROUP  2 

1.  NODES  Operation 

The  NODES  operation  creates  the  matriz  of  nodal 
coordinates  for  a  structural  problen.  This  information  was 
entered  in  cartesian  coordinates,  node  by  node. 

The  user  aay  now  opt  to  enter  data  in  cartesian, 
cylindrical  (any  azis  longitudinal)  or  spherical  coordi¬ 
nates,  It  will  be  converted  to  the  cartesian  coordinate 
system  used  by  CAL.  Additionally,  the  user  aay  generate  new 
nodes  by  specifying  one  node  and  a  node  generation  vector. 
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Y.  SO  LOTI  08  OP  C0BDQCTI08  HEAT  TBABSFBB  PBOBLEHS 

In  order  to  solve  a  conduction  heat  transfer  problem, 
the  user  aust  provide  the  solution  algorithm  to  CAL.  A 
discussion  of  the  matrix  manipulation  and  equation  solving 
techniques  can  be  found  in  Beference  4  .  This  chapter  will 
present  possible  algorithms  for  solution  of  this  class  of 
problems. 

A.  STEADY  STATS  PHOBLEHS 

These  problems  take  tha  fora: 

men  ♦  tF>  =  to j 

1.  Linear  Heat  Conduct  ion  Problems 

This  is  the  siaplest  case  to  consider.  After  input¬ 
ting  mesh  data,  the  user  must  fora  the  conductivity  matrix 
[K],  (SYHC)  and  the  load  vector  {F}  ,  (FORK).  Then  the  nodal 
temperatures  aust  be  calculated  (CALC)  and  printed  (PTEMP). 

Consequently  the  sequence  of  solution  operations  for 
this  type  of  problem  would  be: 

SYHC 
FOB  H 
CALC 
PTE  HP 

2.  Hon-llnear  Eg at  con  duct  ion  Problems 

Since  the  conductivity  matrix  is  time  dependent  on 
temperature,  an  iterative  algorithm  must  be  used.  This 
requires  a  looping  operation  (LOOP,  NEXT)  around  the  linear 
steady  state  sequence. 
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The  solution  operations  say  be: 

LOOP,  II 
STS  C 
FOBS 
CALC 
PTE  BP 
NEXT 

where  N1  is  the  user's  guess  of  the  nuabex  of  iterations 
necessary  to  obtain  equilibrium.  However,  the  prograa  main¬ 
tains  an  internal  check  on  the  residuals.  When  they 
decrease  below  the  predefined  tolerance  (TOL  command, 
default  is  10“*),  the  looping  operation  is  terminated  upon 
the  subsequent  NEXT  command. 

B.  TIHE  DEPENDENT  PROBLEMS 

These  problems  involve  the  full  fora  of  equation  (1)  : 

men  ♦  [c]{T}  ♦  (F)  *  {0}  (i) 

1.  tineas,  flsal  £oMu£t.£on  Eiobleas 

This  case  requires  the  solution  of  a  first  order 
ordinary  differential  equation  (ODE).  Additionally  a  heat 
capacity  natrix,  [ C ] ,  must  be  formed  (CCAP  or  LCAP)  and  a 
time  step  provided  (DTIH) . 

The  differential  equation  solver  is  accessed  using 
the  operation  ODE, HI  where  HI  is  one  of  the  following: 

I NIT  to  specify  initial  temperature  vector  and 
the  integration  constants 

LINE  to  perform  the  two-time  level  algorithm 
QUAD  to  perform  the  three-time  level  algorithm 
The  heat  capacity  and  conductivity  matrices  are 
unchanged  in  a  linear  problem  and  must  be  placed  outside  the 
loop.  The  load  vector  is  reformed  every  time  step  and  the 


These  operations  are 


time  aust  be  advanced  with  ADTia. 
included  in  the  loop. 

The  sequence  of  operations  to  solve  a  linear  tiae 
dependent  problea  aay  be: 

DTI  a 

ODE,INZT 

siac 

CCAP  (or  LCAP) 

LOOP,  SI 
FOB  a 

ODE  ,  LINE  (or  QO  AD) 

ADTIB 
PTE  BP 
NEXT 

where  N1  is  the  nuaber  of  tiae  steps  the  user  wants  to  take. 

2.  Non-linear  Heat  Conduction  Probleas 

The  heat  capacity  and/or  conductivity  matrices  are 
teaperature  dependent  in  this  class  of  problea.  The  temper¬ 
ature  dependent  aatrix  (aatrices)  aust  appear  inside  the 
looping  operation  whereas  the  constant  property  matrix  would 
be  excluded  from  the  loop. 

To  solve  a  fully  non-linear  problem  the  following 
operation  sequence  aay  be  used: 

DTI  H 

ODE, I NIT 
LOOP,  N1 
siac 

CCA P  (or  LCAP) 

FOB  H 

ODE, LINE  (or  QUAD) 

ADTIM 
PTE  HP 
NEXT 

where  N1  is  the  number  of  time  steps  the  user  wants  to  take. 
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C.  IOHEBICAL  EXAMPLES 


1.  Hollow  Cylinder  wi£h  Circumferential  Heating  Strips 

1  hollow  cylinder  with  a  four  inch  outer  diameter 
and  a  three  inch  inner  diameter  was  subjected  to  an  axial 
forced  connvection  condition  in  a  wind  tunnel  by  Professor 
P.  P.  Pucci  of  the  Haval  Postgraduate  School.  There  were 
sixteen  one-quarter  inch  wide  heating  strips  equally  spaced 
ower  180  degrees  of  the  outer  surface  as  illustrated  in 
Pig  are  1. 

The  heating  strips  were  maintained  at  a  constant 
temperature  of  160°  F  and  measurements  of  the  surface  temp¬ 
eratures  between  the  strips  were  made  using  teledeltos 
paper.  The  ambient  temperature  was  60°  F.  The  cylinder  was 
considered  to  have  a  density  of  70  lb/ft3  and  a  specific 
heat  of  0.6  BTO/lbmOF.  Tests  were  made  with  four  heat 
transfer  coefficients  (h). 

The  finite  element  model  took  advantage  of  the 
cylinders  symmetry,  consisting  of  one  half  of  a  heating 
strip  plus  one  half  of  the  interval  between  strips.  The 

element  mesh, as  generated  by  the  PLHX  operation,  is  shown  in 

Figure  2.  The  comparison  of  the  model  data  to  the  tele¬ 
deltos  paper  measurements  was  very  favorable  and  is  shown  in 
Figure  3. 

2.  Transient  Sur face  Temperatures  in  an  Infinite  Plate 

in  infinite  flat  plate  was  considered  as  a  test 

problem  (Figure  4)  .  The  plate  was  at  an  initial  uniform 

temperature  greater  than  the  ambient  temperature  and  then 
exposed  to  convection  conditions  with  a  constant  external 
heat  transfer  coefficient  (h). 

The  approximate  solution  was  obtained  using  a 
Heisler  chart.  The  temperature  of  the  outer  wall  was 
computed  at  1  minute,  5  minutes  and  every  5  minutes  there 


after  up  to  30  Minutes.  The  comparison  between  these  values 
and  those  generated  by  CM  was  very  close  and  is  shown  in 
Figure  5. 
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71.  COHCLOSIOH  S  18  D  HEO  BHEIDATI08S 


The  code  that  was  integrated  into  CAL  provides  an  accu¬ 
rate  and  reliable  Beans  for  solving  a  variety  of  conduction 
heat  transfer  problems.  The  systea  is  user  friendly  both  in 
prompting  for  input  and  detecting  errors.  The  use  of  the 
heat  transfer  group  of  commands  is  encouraged,  as  well  as 
efforts  to  increase  its  versatility. 

Shile  the  capabilities  of  the  prograa  are  significant, 
there  is  rooa  for  iaproveaent.  The  present  version  uses 
priaarily  an  in-core  solution  technique,  which  restricts  the 
problea  size  to  within  the  user's  virtual  aachine  space. 
The  capacity  to  handle  larger  probleas  nay  be  increased 
through  the  use  of  an  out  of  core  technique  for  building  and 
storing  the  conductance  and  capaitance  aatrices,  as  well  as 
equation  solving. 

The  graphics  package  could  be  expanded  and  iaproved  to 
provide  acre  information  to  the  user.  The  plotting  of 
isoparametric  elements  needs  to  include  the  capability  to 
generate  and  plot  curved  lines.  mother  desirable  capa¬ 
bility  is  to  portray  all  three  axes  on  the  screen  with 
apparent  depth  and  allow  rotation  to  any  desired  view.  A 
capabilty  to  plot  isotherms  could  be  added  to  enhance  the 
interpretation  of  the  results. 
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Figure  2.  Graphic  Display 
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Infinite  Slab 
Figure  4. 
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APPENDIX  A 


USER'S  GUIDE 

This  appendix  provides  details  on  the  use  of  CAL  with 
the  IBH  3033  computer  at  NPS.  The  prograa  as  nodified  at 
the  Naval  Postgraduate  School  is  subseguently  referred  to  as 
CAL-NPS.  Section  A  provides  the  details  on  the  coaaand 
structure.  Section  B  is  a  summary  of  coaaands  available. 
Section  C  provides  the  job  control  language  for  executing 
the  prograa  in  both  the  batch  and  interactive  aodes  at  NPS. 
Section  D  contains  detailed  specifications  for  each  avail¬ 
able  command.  Finally,  section  E  gives  direction  for 
solving  larger  probleas  with  CAL-NPS.  The  majority  of  this 
appendix  was  originally  published  as  Reference  1  .  The 
author  wishes  to  express  appreciation  to  Professor  Wilson 
for  peraission  to  use  this  material. 

A.  FORM  AID  RESTRICTION  OF  THE  LANGUAGE 

CAL-NPS  is  an  interpretive  language  which  is  designed  to 
manipulate  arrays  and  matrices,  to  perform  standard  struc¬ 
tural  analysis  operations  and  to  perform  conduction  heat 
transfer  analysis  operations.  A  CAL-NPS  program  run 
involves  the  reading  of  the  input  deck  once  and  executing 
the  commands  designated  by  the  operation  cards  as  they  are 
encountered.  Looping  operations  allow  a  sequence  of 
commands  to  be  executed  more  than  once. 

The  input  deck  is  composed  of  operation  cards  and  data 
cards.  The  data  cards  directly  follow  each  operation  card 
which  requires  data  (see  LOOP  operation  for  exception  to 
this) .  The  operation  card  contains  the  name  of  the  opera¬ 
tion  to  be  executed,  names  of  arrays  associated  with  the 
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operation  and  integer  constants.  Examples  of  the  general 
fora  of  this  card  are: 

OP,H1,H2,H3,H4,H5  ,H1,  N2 , M3 ,  H4  comments 
OP, HI, HI, H2 
OP,  HI 
OP 

in  which  OP  is  the  name  of  the  operation  to  be  executed.  Hi 
is  the  name  of  an  array  and  Hi  is  an  integer.  The  names  of 
OP  or  Hi  are  one  to  eight  alphabetic  or  numeric  characters 
to  be  selected  by  the  user.  The  first  character  of  a  name 
must  be  alphabetic.  The  sequence  of  terms  OP,  Hi  and  Hi 
must  be  separated  by  commas.  Characters  following  a  blank 
will  be  printed  as  comments  in  the  output  from  the  program 
run. 

If  an  operation  attempts  to  load  or  generate  an  array 
which  previously  existed,  the  program  will  delete  the  array 
before  the  execution  of  the  operation,  k  new  array  need  not 
be  the  same  size  of  the  old  array  which  had  the  same  name. 
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B.  SOHHARY  OP  COHM1NDS 


1.  SaaSIfii  gaaj&Ms 


*  indicates  a  significant  change  or  addtion  in  CAL-NPS 

START  -  Initialize  for  the  next  problem 

STOP  -  Normal  termination 

NO  -  Temporary  suppression  of  output 

YES  -  Restores  output 

LABEL  -  Print  comments 

READ  -  Change  logical  device  for  input 

BRITE  -  Change  logical  device  for  output 

TIBS  -  Suppress  tine  printout 

SAVE  -  Interrupt  a  problem 

RESOHE  -  Continue  an  interrupted  problem 

LIST  -  List  arrays  and  storage  used 

HELP  -  Access  HELP  files 

2.  General  Matrix  Commands 

LOAD  -  Load  user  defined  real  matrix 

LOADI  -  Load  user  defined  integer  matrix 

ZERO  -  Create  null  or  unit  matrix 

PRINT  -  Matrix  print  operation 

DUP  -  Matrix  duplication 

ADD  -  Matrix  addition 

SOB  -  Matrix  subtraction 

MOLT  -  Matrix  multiplication 

TRAN  -  Matrix  transpose 

SCALE  -  Multiply  a  matrix  by  a  scalar 

SOLVE  -  Solution  of  linear  equations 

DOPSH  -  Form  sub-matrix  from  large  matrix 

STOSM  -  Store  sub-matrix  in  large  matrix 

DOPDG  -  Form  row  matrix  from  diagonal 

STODG  -  Store  row  on  diagonal 
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BiX  -  Evaluate  row  aaxiauas 

HOB H  -  Evaluate  aatrix  norms 

INTEL  -  Invert  each  term  in  aatrix 

SQREL  -  Square  root  of  each  term  in  aatrix 

LOG  -  natural  log  of  each  tera  in  matrix 

PROD  -  Evaluate  product  of  all  terms  in  a  aatrix 

DELETE  -  Delete  aatrix  fro  a  storage 

3.  Static  Analysis  Operations 

NODES  -  Input  structural  joint  geometry 

BOUND  -  Specify  boundary  conditions 

BEAN  -  Fora  3-D  beam  stiffness  aatrix 

TRUSS  -  Fora  3-D  truss  stiffness  aatrix 

PLANE  -  Fora  3  to  8  node  plane  stiffness  aatrix 

SLOPE  -  Fora  stiffness  aatrix  froa  slope/deflection  eq. 

FRA  HE  -  Fora  2-D  fraae  stiffness  aatrix 

LOADS  -  Fora  load  vector 

ADDSF  -  Fora  global  stiffness  and  mass  matrices 

ADDK  -  Add  element  aatrix  to  global  aatrix 

NEHFRC  -  Calculate  element  forces  from  joint  displacements 

DI5PL  -  Print  joint  displacements 

FORCE  -  Evaluate  and  print  member  forces 

fixuais  s  SfigsaUsBs 

FUNG  -  Generate  equal  interval  time  function 
STEP  -  Integrate  dynamic  equilibrium  equations 
EIGEN  -  Evaluate  mode  shapes  and  frequencies 
DINAH  -  Evaluate  uncoupled  equations  of  motion  by 
mode  superposition  nethod 
PLOT  -  Line  printer  plot  of  joint  time  history 

5.  geat  Transfer  Operations 

HTXFR  -  Initiate  heat  transfer  problem 

COORD  -  Input  nodal  coordinates 

ELCON  -  Input  element  connectivity  aatrix 
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PROP  -  Input  aaterial  property  data  * 

CTBHP  -  Input  constant  teaperature  node  data  * 

PROP  -  Bstablish  profile  of  equations  * 

SYHC  -  Create  syaaetric  conductance  aatrix  * 

OSIHC  -  Create  unsyaaetric  conductance  aatrix  * 

LCAP  -  Create  luaped  capacitance  aatrix  * 

CCAP  -  Create  consistent  capacitance  aatrix  * 

FORM  -  Create  flux  rector  * 

CALC  -  Solve  tiae  independent  systeas  of  equations  * 

ODE  -  solve  tiae  dependent  systeas  of  equations  * 

PTE HP  -  Print  nodal  teaperatures  * 

TOL  -  Set  solution  convergence  tolerance  * 

CON?  -  Perfora  teaperature  convergence  test  ■* 

DTIH  -  Set  tiae  step  increaent  * 

ADTIH  -  Advance  tiae  by  one  tiae  step  * 

PROMPT  -  Suppress/restore  proapts  * 

6 .  Graphics  Operations 

GRAPH  -  Initiate  graphics  * 

TITLE  -  Label  aesh  plot  a 

PLH X  -  Plot  heat  transfer  aesh  * 

PLST  -  Plot  structural  aesh  * 

7.  Loop  Operations 

LOOP  -  Start  of  loop 

HEXT  -  End  of  loop 

SKIP  -  Conditional  skip  of  operations  within  loop 

8.  Bases  Available  fat  Peer  Spbcputlnes 

OSERA 

□SERB 
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C«  EXECUTION  OF  CIL-NPS 


Por  the  tiae  sharing  (CHS)  system  at  Naval  Postgraduate 
School,  do  the  following: 

_(0se  the  standard  LOGON  procedue) 

_link  0040P  191  199 
ENTER  PASS WORD : 

_XXXXX  (ESAN) 

R;  T=0.01/0.01  11:09:55 

^access  199  C  (Note:  You  must  access  the  nC" 

C  (199)  H/O  disk) 

B;  T®0.01/0. 01  11:10:04 
_cal 

ENTER  TE RHINAL  CODE: 

1  *  PLOT- 10  Compatible  Terminal  (GRAPHICS) 

2  *  IBH  3277  DUAL  SCREEN  (GRAPHICS) 

3  =  Any  Alpha  Numeric  Terminal  (NO  GRAPHICS) 

_1  (or  2  or  3,  as  appropriate) 

(The  computer  responds  with  several  lines  of  procedure) 
EXECUTION  BEGINS  . 

- 0.0  SECONDS 

(You  are  now  under  the  control  of  CAL-NPS) 

start 

** START 

) 

(Your  own  CAL-NPS  program  is  inserted  here) 

stop 

eesTOP 

R;  T*0. 01/0. 0 1  11:12:45 
_log  (Terminates  session) 
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D.  CAL -UPS  COHHAHD  SPECIFICATIONS 


1.  General  Matrix  Operations 

CAL-MPS  has  most  of  the  standard  matrix  operations 
plus  some  special  array  operations  which  are  useful  in  engi¬ 
neering  analysis.  The  following  is  a  list  of  approximately 
32  operations  which  are  used  for  control  and  general  aatrix 
■an ip alation. 

A  indicates  the  foraation  of  a  new  aatrix.  A 
aatrix  previously  defined  with  the  same  name  will  be 
deleted.  A  indicates  modification  of  an  existing 

matrix. 

Mote:  Whenever  the  expression  ’’card"  is  used  it  is 

meant  to  also  stand  for  ’’instruction”  in  interactive  mode. 


START 


This  operation  eliainat 
loaded  or  generated. 


es 


all  arrays 


which 


were  previously 


STOP 

This  operation  causes  normal  termination  of  a  CAL-NPS 
program. 


NO 

YES 


These  operations  are  used  to  selectively  suppress  output 
from  CAL-MPS.  The  MO  operation  suppresses  all  printing, 
except  diagnostics,  until  the  operation  YES  is  encountered. 
Therefore,  in  subsequent  runs  or  the  same  CAL-MPS  program, 
output  which  was  previously  correct  need  not  be  reprinted  it 
these  cards  are  inserted  in  the  data  deck. 


LABEL.N1 


.s  operation  wj.ll  read  and  print  Ml  comment  cards,  which 

_ .low  the  operation  card.  Column  1  of  each  card  will  be 

interpreted  as  a  standard  carriage  control  symbol  (i.e.  0 

for  double  space  and  1  for  skip  to  the  top  of  the  next 
page)  . 
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HEAD, HI 

THIS  QPEHATIOH  IS  VALID  OH LY  WITH  THE  CP/CMS  TIME  SHAEIHG 
SI?  TEH. 

This  operation  permits  the  selection  of  a  user's  file  or  the 
terainal  as  the  input  file  device.  The  default  is  the 
terainal.  If  HI  is  4,  subsequent  commands  will  be  read  froa 
FILE  FT04F001  on  the  user's  A-disk.  If  HI  is  5,  the 

terainal  will  be  restored  as  the  input  file  device.  All 
dish  fi4.es  prepared  for  use  with  this  command  should  end 
with  either  STOP  or  EEAD,5.  This  command  will  not  be 
executed  on  the  MVS  (batch)  system. 


KBITS* HI 

THIS  OPERATION  IS  VALID  OH  LY  KITH  THE  CP/CHS  TIME  SHAEIHG 
SY?TEH. 

This  operation  permits  the  selection  of  a  file  on  the  user's 
A-disk  or  the  terainal  as  the  output  file  device.  The 
default  is  the  terminal  and  all  error  message?  will  be 
printed  at  the  terminal  regardless  of  the  output  file  device 
selected.  If  Hi  is  8,  subsequent  ouput  will  be  written  into 
FILE  FT08F00 1  on  the  user's  A-disk.  If  HI  is  6,  the 
terminal  will  be  restored  as  the  output  file  device.  This 
coaaand  will  not  be  executed  on  the  MVS  (batch)  system. 


TIME 

This  operation  permits  the  time  printout  to  be  suppressed 
without  los?  of  other  output.  A  second  TIME  will  restore 
the  time  printout  unless  the  print  output  is  suppressed  with 
the  HO  coaaand. 


SAVE  or  SAVE, HI  or  SA VE, Ml , HI 

This  operation  saves  all  arrays  in  storage  at  the  time  of 
issuance.  Saved  arrays  will  contain  all  modification?  made 
since  their  creation.  Ml  and  HI  are  optional  and  if  not 
included,  the  arrays  will  be  stored  in  FILE  02  on  the  user's 
A-disk.  The  saved  files  will  be  assumed  to  be  for  general 
matrix  manipulation  or  a  structural  problem.  Ml  is  the  file 
naae  (up  to  six  letters)  under  which  the  user  wishes  to 
store  the  arrays.  The  file  type  will  be  SAVE.  If  Hi  is  1, 
a  general  matrix  manipulation  or  structural  problem  is  being 
saved.  If  HI  is  2,  a  heat  transfer  problem  is  being  saved. 


RESUME  or  RESUME, HI  or  EESUME,M1,N1 

This  operation  reads  a  saved  file  into  memory..  Any  arrays 
currently  in  storage  will  be  destroyed.  A  file  must  have 
been  previously  created  on  the  user's  A-disk  using  the  SAVE 
operation.  Ml  and  Hi  are  optional  and  if  not  included,  FILE 
02  will  be  read.  It  will  be  assumed  to  be  a  general  matrix 
manipulation  or  structural  problem.  Ml  is  the  file  name 
assigned  to  the  saved  file  on  the  user's  A-disk.  The  file 
type  must  be  SAVE.  If  HI  is  1 ,  a  general  matrix  manipula¬ 
tion  ox  structural  problem  is  being  resumed.  If  Hi  is  2,  a 
heat  transfer  problem  is  being  resumed. 


LIST  .  ... 

The  LIST  operation  prints  the  directorv  information  for 
arrays  in  storage  and  the  amount  of  storage  used. 
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LOAD,H1,N1,H2,H3 

This  operation  will  load  a 
which  has  HI  rows  ana  N2  c 
ar$  punched  in  row-wj.se  s 
this  operation.  N3  is  opt 
the  cards  are  punched  in  (8? 
additional  card  which  conta 
oast  proceed  the  data,  Fo 
numbers  per  card  in  field  w 
would  contain  the  following 
nine,  the  data  cards  will  5 


n  array  of  real  nuabers  named  HI 
olumns,  The  teras  of  the  array 
gguence  on  data  cards  following 
lonal.  If  N3  is  zero  or  blank, 
10.0)  format.  If  N3  is  one,  an 
ins  the  format  of  the  data  cards 
r  example,  if  the  data  is  to  be  4 
ldths  of  15,  the  additional  card 
information:  (4P15.0)  .  If  H3  is 
e  read  in  free  format. 


LOA DI,H1,N1,N2,N3,N4  or 
LOA  01, HI, HI, N2, M3  or 

LOADI,H1 ,H2,N1, N2,N3, N4 


L0ADI,H1  ,111  ,H2,  ,H4 
LOA  DI,H1 , N1 ,N2  or 


or 


integer  array  named  Hi  which 


IU  AUUVMVfc  i d V 

The  terms ^of. the  array  are  punc 


611 


This  operation  will  load  an 

91  rows  ana  N2  columns.  T _ _  __  _ _ _  _ c _ T _ 

in  row  wise  sequence  on  data  cards  following  this  operation. 
H2, H3  and  N4  are  optional.  If  H3  is  zero  or  blank,  the  data 
must  be  punched  in  (1615)  format.  If  N3  is  one,  an  addi¬ 
tional  card  containing  the  format  of  the  data  cards  must 
follow  this  operation  and  preceed  the  data.  For  example,  if 
the  data  is  to  be  4  numbers  per  card  in  field  widths  of  10, 
the  additional  card  would  contain:  (4110).  If  N3  is  nine, 
the  data  will  be  read  in  free  format. 

If  the  letters  "INCH"  are  placed  m  the  position  of  H2,  this 
operation  has  an  increment  generation  capability.  Data  must 
be  entered  as  follows: 

Item  Contents 

1  Row  number 

2  Value  1 

3  Value  2 


H2+1  ‘Value  N2 

H2+2  Generation  code 

If  the  generation  code  is  not  zero,  the  next  card  must 
contain  the  following: 

Item  Contents 

1  Row  number  increment 

2  Value  1  increment 

3  Value  2  increment 

N2+1  ‘value  N2  increment 
N2+2  Last  row  to  be  generated 


ZERO,  H1,N1,N2,N3,N4 


A  real  matrix  named.  HI  is 
columns.  The  terms  m  this 
values: 


Therefore 

matrices. 


HI 

HI 

th 


11:51 

is  c 


13 

N4 


I  = 
J  = 


operation  can 


and  . N2 
ollowmg 


1 ,  .  .  .  Ml 
1.. ..N2 

be  used  to  form  null  or  unit 
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PRINT, HI  or  PRINT, HI, HI  or  PRINT, HI  r  HI,  M2  or 

PRINT,  HI,  Nt,N2,N3 


This  operation  will  print 
format  of  up  to  eight  colun 
are  optional.  N1  is  thg  n 
the  operation  card)  which 
defaults  to  zero.  If  N2 
printed  in  partitioned  for 
Lines  will  have  N2*15  ♦  5 
printing  125  charcters  per 
greater  than  zero,  integer 
value  is  zero  and  real  fora 
cautioned  not  to  overcone 
device  in  use  to  avoid  wrap 


the  array  naaed  HI  in  a  aatrix 
ns  per  line.  H1,N2,N3.  and  NQ 
umber  of  consent  cards  (following 
will  be  read  and  printed.  Hi 
is  included,  the  aatrix  will  be 
a  with  N2  coluans  per  partition. 

characters.  N2  defaults  to  8, 
line.  If  N3  is  included  and  is 
foraat  (16)  is  used.  The  default 
at  (PD15.7)  is  used.  The  user  is 
the  capacity  of  the  displaying 
around  on  trie  screen. 


♦ 

DUP  ,  Ml,  M2 

This,  operation  will  fora  an  array  naaed  H2  which  is  iden¬ 
tical  to  the  array  naaed  HI . 


ADD  ,H 1, H2 

This  operation  will  replace  aatrix  HI  with  the  sun  of  the 
eat  rices  HI  and  H2. 


SUB  ,H1,  H2 

This  operation  will  replace  matrix  HI  with  aatrix  Hi  less 
matrix  H 2. 


HULT, fll , H2 , H3 

This  operation,  generates  a  new  aatrix  H3  which  is  the 
product  of  matrices  Hi  and  H2,  or  H3  =  H1ftH2. 


♦ 

TRAN,  HI,  H2 

This  operation  generates  a  new  matrix  H2  which  is  the  trans¬ 
pose  of  matrix  HI. 


_  # 
SCALE, HI, H2 


This  operation  replaces  each  term  in 
with  the  term  multiplied  by  the  term 
naaed  H2. 


the  natr 
H2  (1,1)  0 


i 


x  named  Hi 
the  matrix 
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SOLVE , Ml ,H2,N1,M2  or  SOLVE,  Ml ,M 2, ,  N2  or  SOLVE, Ml , HI , M2  or 
SOLVE, Ml, M2, HI  or  SOLVE, HI, M2 


If  HI  a  0#  this  operation  solves  the  matrix  equation  AX=B. 
Hi  is  the  name  of  the  &  natrix  ana  M2  is  the  name  of  the  B 
matrix.  Hatnx  A  is  triangularized  and  the  results,  X,  are 
stored  in  M2.  If  111  =  1,  Matrix  &  is  triangularized  only. 


stored  m  M2,  if  ni  =  i,  Matrix  k  is  triangularized  only. 

If, M1*2,  for  a  given  B  aatrix  and  the  A  matrix  previously 
triangularized.  tne  B  matrix  is  replaced  by  the  results,  X. 
If  M  1*3,  Matrix  k  is  replaced  by  its  inverse  POE  SYMMETRIC 
MATRICES  ONLY. 

If  H2*0  or  blank,  matrix  k  is  symaetric.  If  M2  is  nonzero 
the  aatrix  A  is  not  symaetric. 

For  symmetric  matrices,  A  is  factored  into  the  LDL  fora. 
The  diagonal  D  matrix  is  stored  on  the  diagonal  of  A.  The 
parameter  M2  permits  the  direct  solution  of  non-syametric 
systems  of  equations.  If  M2  is  not  equal  to  zero,  an  LO 
decomposition  of  matrix  A  will  be  performed.  Mo  direct 
replacement  of  Ml  by  its  inverse  is  available  for  the  non- 
sym metric  case.  Instead,  use  the  ZERO  operation  to  create 
an  identity  aatrix  M2  of  the  same  order  as  HI.  The  command 
SOLVE, Ml , M2,, M2  will  then  replace  the  aatrix  M2  with  the 
inverse  of  the  aatrix  A. 


DUPSM ,M 1 , M2, Ml,  M2,  M3,  M4 


This  operation  forms  a  new  subaatrij:  named  M2  with  M3  rows 
ana  N4  columns  from  the  terms  within  the  matrix  named  Ml. 
The  first  term  of  aatrix  M2,  112(1.1),  will  be  from  row  N1 
and  column  N2  of  aatrix  Ml,  or  Ml  (Ml,  M2). 


ST0SM,H1,M2,N1,N2 

This  operation  stores  a  subaatrix  named  M2  within  the  matrix 
named  Hi.  The  first  term  of  the  subaatrix  M2  will  be  stored 
at  row  Ml  and  column  M2  of  matrix  Ml.  The  terms  within  the 
area  of  Ml  in  which  M2  is  stored  will  be  destroyed. 


D0PDG,M1,M2 

This  operation  forms  a  new  natrix  named  M2  from  the  diagonal 
terms  of  Ml. 


STODG  ,M1  ,M2 

This  operation  stores  a  row  or  column  matrix  named  m2  at  the 
diagonal  locations  of  natrix  Ml. 


MAX  , M  1,  M2 


This  operation  forms  a  column  natrix  named  M2 
row  contains  the  maxiaua  absolute  value  of  the 
row  in  matrix  Ml.  The  maximum  and  its  col 
printed  for  each  row. 


M2  in  which  each 
the  corresponding 
column  number  is 


HORN, Ml, M2, HI 

If  HI  *  0,  a  row  aatrix  naned  M2  is  formed  in  which  each 
column  contains  the  sum  of  the  absolute  values  of  the  corre¬ 
sponding  column  of  matrix  si.  If  N1  is  not  equal  to  0,  a 
row  matrix  named  H2  is  formed  m  which  each  column  contains 
the  square  root  of  the  sum  of  the  squares  of  the  values  of 
the  corresponding  columns  of  aatrix  HI. 


SQREL,H1 

This  operation  replaces  each  tera  in  the  aatrix  named  Ml 
with  the  square  root  of  the  tera. 


LOG, Ml 

This  operation  , replaces  each  tern  in  the  matrix  named  Ml 
with  the  natural  log  or  the  term. 


PROD,  Ml,  M2 

This  operation  forms  a  1  x  ?  array  named  M2  which  contains 
the  product  of  all  terms  m  the  matrix  naned  Ml.  The 
product,  Z,  is  stored  as  two  numbers  of  the  form: 

X  *  P*10*«E 

in  which  M2(1)  *  P  and  M2(2)  =  E,  the  exponent. 


DELETE,  Ml 

This  operation  will  cause  the  elimination  from  storage  of 
the  array  naned  HI. 


CMS, HI 


This  operation  allows  the  user  to  issue  CMS  commands  while 
under  the  control  of  CAL-HPS.  In  general  a  command  that 
reloads  the  virtual  core  will  not  be  allowed.  Examples  are 
FORTH!.  SCRIPT,  XEDIT,  any  language  processor,  SORT,  and 
other  large  system  nodules  such  as  COPYFILE  and  MOVFILE.  N1 
is  the  number  of  words  in  the  command  (1  to  9  are  allowed!. 
Hote  that  a  parenthesis  -  or  ") "  is  counted  as  a  word. 
All  words  must  be  left  justified. 


fhe  user  will  be  prompted  for  each  word.  If  this  operation 
s  used  in  an  PT09F001  FILE,  each  word  nr'-"1-  --  -  - 


lin  e. 


FILE,  each  word  must  be  on  a  separate 
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2.  ££&ti£  Afl&jxsis  2E£taiians 


The  purpose  of  this  series  of  operations  is  to  for» 
the  total  stiffness  and  diagonal  (luaped)  nass  aatrices  for 
systees  of  two-  or  three-dimensional  elements.  Por  three- 
disensional  analysis  there  are  beam  and  truss  elements 
available.  For  two-dimens ional  analysis,  there  is  a  frame 
element,  a  slope/deflection  eleaent  for  beams,  and  a  4  to  8 
node  isoparametric  finite  element  available. 

After  the  creation  of  an  array  containing  the  coor¬ 
dinates  of  the  joints  of  the  system,  the  specification  of 
displacement  boundary  conditions,  the  tabulation  of  material 
and  section  properties,  the  mass  and  stiffness  aatrices  are 
formed  for  each  structural  member  and  placed  in  sequence  on 
low  speed  storage  along  with  the  global  equation  numbers 
which  are  associated  with  their  stiffness  terms.  In  addi¬ 
tion,  the  member  force-displacement  transformation  aatrices 
are  formed  and  stored  on  a  separate  low  speed  storage  file 
along  with  the  appropriate  displacement  numbers. 

The  NODES  operation  is  used  to  specify  or  generate 
the  geometry  of  the  system.  The  operation  BOUND  specifies 
which  joint  displacements  exist  and  assigns  internal  equa¬ 
tion  numbers  to  these  displacements.  Therefore,  each  joint 
may  have  from  zero  to  six  displacement  degrees  of  freedom. 
Tables  of  material  and  section  properties  for  the  various 
members  are  loaded  and  printed  as  standard  arrays  of 
information. 

&  special  operation,  ADDSF,  is  used  for  the  direct 
addition. of  element  stiffnesses  to  form  the  total  stiffness 
and  diagonal  nass  matrix  of  the  system.  The  ADDK  operation 
may  be  used  to  add  individual  elements  into  the  total  system 
aatrices.  The  LOADS  operation  specifies  the  concentrated 
joint  loads  for  all  load  conditions.  After  the  direct  solu¬ 
tion  for  joint  displacements  due  to  static  or  dynamic  loads. 
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the  eeeber  forces  can  be  evaluated  using  the  (iEMFBC  opera- 
tion.  The  DISPL  operation  is  used  to  print  the  displace- 
aents  in  joint  number  order. 
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♦  ♦ 

NODES ,H 1  ,N  1  or  NODES, Ml, N1, N2 


This  operation  generates  a 
contains  tne  coordinates 
system.  N2  is  optional, 
for sat  as  folio vs; 


lor1!!!*  folnl’s 
Da±a  must  be 


named  Ml  which 
m  a  structural 
entered  in  free 


Item  Contents 

1  Node  number 

2  X-coordinate 

3  Y-  coordinate 

4  Z-coordmate 


If  N2  *  1,  there  i$  a  joint  generation  and  coordinate  system 
conversion  capability.  Data  must  be  entered  in  free  format 
as  follows; 


Item  Contents 

1  Node  number 

2  X-coordinate 

3  Y- coordinate 

4  2-coordinate 

5  System  type 

6  Generation  code 


System  type  refers  to  th?  system  used  when  input ing  the 
data.  All  coordinates  will  be  converted  to  the  cartesian 
system  for  use  by  CAL-NPS. 


System  type 

2 

3 

4 

5 


S  yst  em . 

Cartesian 

Cylindrical,  Z  axis  longitudinal 
Cylindrical,  Y  axis  longitudinal 
Cylindrical,  X  axis  longitudinal 
Spherical 


The  input  data  is  the  same  as  above  with  the  following 
correspondence; 


£2£tesj£n 

Y 

Z 


Cylindrical 

r 

e 

z 


Spherical 

(3 

* 


If  the  generation  code  is  not  zero,  the  next  card  is. a 
generation  vector  for  the  self  generation  of  nodes.  It  is 
formatted  as  follows; 


Item  Contents 

1  Node  number  increment 

2  X  increment 

3  Y  increment 

4  2  increment 

5  Last  node  number  to  be  generated. 

It  is  assumed  that  the  increments  pertain  to  the  same  system 
of  reference  as  the  preceedipg  card. 

This  operation  must  be  terminated  by  a  line  of  alternating 
zeros  and  blanks. 
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BOUND, HI 

?his  operation  specifies  the  displacements  which  are  nonzero 
or  the  structural  system  of  joints  specified  by  the  NODES 
operation.  Where: 

HI  =  Name  of  boundary  condition  code  array  to  be 
generated. 

This  operation. is  followed  by  a  series  of  cards  containing 
the  following  information  m  free  format: 

Item  Contents 

1  Nod$  number  for  the  first  node  in  a 
series  of  nodes  with  identical 
displacement  specification. 

2  Node  number  for  the  last  node  in  the 

series. 

3  X-trans  lation 

4  Y-trans lation 

5  Z-translation 

6  X -rot at  ion 

7  Y-rotation 

8  Z-rotation 

9  Node,  number  i  ncrement  used  to  generate 

conditions  for  addtional  nodes. 

J  translation  or  rotation  equals:  (a)  zero  fon  zero  or  unde- 
med  displacements,  or  (b)  one  for  nonzero  displacements  to 
be  evaluated  by  other  operations. 

If  a  node  boundary  condition  is  not  specified,  all  displace¬ 
ments  at  that  node  are  assumed  zero.  Cards  may  be  supplied 
in  any  order.  If  node  boundary  conditions  are  specified 
more  than  once,  the  last  definition  is  used.  This  sequence 
of  data  must  be  terminated  by  a  card  of  alternating  zeros 
and  blanks. 

She  selection  by  the  user  of  which  nodes  have  nonz?£o 
isplacments  requires  an  understanding  of  the  direct  stiff¬ 
ness  procedure.  Displacement  degrees  of  freedom  which  have 
no  stiffness  associated  with  the  displacement  must  be 
considered  to  be  undefined  since  it  is  not  possible  to 
develop  an  equilibrium  equation  for  that  direction.  The 
total  number  of  nonzerp  displacements  specified  will  be  the 
size  of  the  total  stiffness  matrix  to  be  defined  by  the 
ADDS?  operation. 
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♦ 

BEAM,  Ml  ,M2,M3,M4 


This  operation  calculates  the  element  spiffness,  mass  and 
force-displacement  transformation  matrices  for  3-D  beam 
members.  These  arrays  are  stored  in  sequence  on  low  speed 
storage  to  be  used  by  other  operations  where: 


HI  is  the  name  of  the 
M2  is  the  name  of  the 


_  _  le  beam  element  group 

is  the  name  of  the  coordinate  array 
M3  is  the  name  of  the  boundary  condition  array 
M4  is  the  name  of  the  array  which  contains  beam 
properties  and  has  been  loaded  by  the  standard 
matrix  LOAD  operation 


Qne  card  for  each,  beam  in.  this  group  of  beam  elements  must 
tollow  this  operation.  The  Beam  cards  are  punched  m  free 
format,  where: 

Item  Contents 

1  Beam  identification  number 

2  Node  number  1 

3  Node  number  <1 

„  4  Node  number  K 

5  Beam  property  number  NP 

This  sequence  of  cards  must  be  terminated  with  a  card  of 
alternating  zeros  ana  blanks. 


The  material  and  geometric  properties  for  each  element  are 
given  in  the  H4  array  in  the  following  order: 


N4 

M4 

H4 

MU 

H4 

H4 

M4 


NP  ,  1 
NP,  2 
NP  ,  3{ 
NP,4‘ 
NP,  5 
NP,  6 
NP  *  7J 


Axial  area  of  member,  A 

Torsional  moment  of  intertia,  J 

Moment  of  inertia  about  axis  2,  I 

Moment  of  inertia  about  axis  3,  I 

Modulus  of  elasticity,  E 

Shear  modulus,  G 

Mass  per  unit  length  of  beam 


phere  NP  is  the  specific  material  propert- 
m  item  5  of  the  beam  c$rd .  The  local 
given  in  the  following  figure. 


. number  specified 
ign  convention  is 
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DEFINITION  OF  POSITIVE  BEAM  FORCES 
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TBOSS'BI  ,82, S3,  H4 


■ass  and  force- 
trass  Benners, 
ia  sequence  and 


the  aaae  of  this  group 
the  naae  of  the  cdardii 


of  trass  aenbers 

-ir- _ coordinate  array 

s  the  naae  of  the  boundary. condition  array 
_.an  IP  by  3  array  or  section  properties  la 
which  HP  is  the  naaber  of  different  section 
properties  and 
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84 

84 


(HP  ,  1)  * 
Hp;2  » 
HP, 3  - 


The  cross-sectional  area,  \ 

The  aodalas  of  elasticity,  E 
the  aass  per  unit  length  of  the  aeaber. 


This  aatrix  can  be  loaded  by  the  aatrix  LOAD  operation. 


This  operation  is  t  followed  by  one  card,  per  truss  aeaber  in 
free  foraat  with  the  following  mioraation: 

I  tea  Contents 

1  Truss  aeaber  identification  naaber 

2  Joint  naaber  I 

3  Joint  naaber  J 

4  Section  property  nuaber,  HP 


This  operation  aast  be  terainated  by  a  card  of  alternating 
zeros  ana  blanks. 


maE, 81,82 


This  operation 
two-diaensional 


for  as  theL  6  x.  6  stiffness 
raae  aeaber  shown  below. 


aatrix  for 


the 


FRAME  MEMBER 


51 


The  properties  of  the  aeaber  are  defined  on  one  cards  iaae- 
diately  following  the  PR1BB  operation  card.  This  second 
card  is  punched  in  free  foraat  and  contains  the  following 
inf  oraation: 


ixiaf  area,  k  a  , 
Hodulus  of  elasticity, 
Sgaent  of  inertia,  I 


B 


GEOMETRY  AND  JOINT  DISPLACEMENTS 


12  is  a  3  i  6  f orce-disp lac eaent  tpansf oraation  aatrix  which 
is  based  on  the  positive  definition  of  the  eleaent  forces 
shown  below. 


These  forces  car  be  calculated  fros 
equation  with  the  SESPBC  operation. 

“u«n 

|U,» 

Um 


the  folio sing 


M 

M 

P 


=  [M2] 


U* 

Us 

U* 


aatriz 


SLOPE, HI 


This  operation  fores  a  4  z  4, stiffness  aatriz.  Hi 
or  coluan  aeeber  tfoa  the  classical  slope-def lection  equa¬ 
tions,  The  properties  of  the  aesber  are  defined  on  one  card 
iaaediately  following  the  operation.  This  second  card  is 
punched  in  free  foraat  and  contains  the  following  inforaa- 
tion: 


fpr  a  beaa 


Item 

1 

2 

3 


Contents 

Bonent  of  inertia,  I 
Bodulus  of  elasticity,  2 
Length  of  aeaber,  L 


The  sign  contention  is  defined  as  follows: 


I—  -*4 

DISPLACEMENTS 


The  aeaber  forces  are  defined  in  t§»ras  of  joint  displace¬ 
ments  by  the  following  slope  deflection  equations. 


V,=  -v,= 
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PLANE, Ml, M2, M3,  H4,N1,  N2 

This  operation  calculates  the.  element  stiffness, lass  and 
stress-displacement  transformation  matrices  for  4  to  8  node 
isoparametric  elements  (Y-Z  plane  only).  These  arrays  are 
stored  <?n  low  speed  storage  to  be  used  later  by  other  opera¬ 
tions  (x.e. ,  ADDSF  and  FORCE).  The  arguments  are  defined  as 

HI  is  the  user  defined  n^ae  of  th?  eleaent  group 
M2  is  the  naae  of  the  joint  coordinate  array 
M3  is  the  naae  of  the  Boundary  condition  array 
M4  is  the  naae  of  the  array  which  contains  the 

saterial  properties  of  the  eleaents  (one  row  per 
different  material)  where 


(NP,  1)  ~  Modulus  of  elasticity,  E 
NP,2  =  Poissons  ratio.  v 
(NP,  3  *  Thickness  of  element 

(NP , 4  *  Mass  density  of  the  element. 

NP  is  the  material  identification  number. 


N 1  ?md  N?  are  the  nupber  of  integration  points  in  the  r  and 
s  directions  respectively. 

?ne  card  for  each  3  to  8  node  element  in  the  group  must 
oilow  the  operation  card.  The  cards  are  punched  it  free 
format  and  contain  the  following  information: 

Item  Contents.  .  . 

1  Element  identification  number 

2  Node  number  N1 

3  Node  number  N2 

4  Node  number  N3 

5  Node  number  N4 

6  Node  number  N5 

7  Node  number  N6 

8  Node  number  N7 

9  Node  number  N8 

10  Material  identification  number,  NP 

11  Natural  Coordinate  of  stress  output  rl 

12  Natural  Coordinate  of  stress  output  si 

13  Natural  Coordinate  of  stress  output  r2 

14  Natural  Coordinate  of  stress  output  s2 

15  Natural  Coordinate  of  stress  output  r3 

16  Natural  Coordinate  of  stress  output  s3 

N4  through  N8  are  optional,  but  zeros  must  be  inserted  for 
them  if  unused.  The  midside  nodes,  if  present,  must  be 
within  the  center  half  of  the  side.  The  local  numbering 
system  for  the  eleaent  is  shown  m  the  following  figure. 
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vai 


ISOPARAMETRIC  element 


arm***' * 


,th*v,£»S h  il!rati?h8aforclsthare 


♦ 

LOADS, Ml, M2, SM 


This  operation  forms  a  load  matrix  ,  naaed  Ml  of  ,  HI,  columns 
(Ni  load  conditions)  where  M2  is  the  naae  or  the  boundary 
condition  array  generated  by  the  operation  BOUND.  This 
operation  is  followed  by  a  series  pr  cards  -  one  for  each 
loaded  joint  for  each  load  condition.  These  cards  are 
punched  m  free  foraat  as  follows: 


Item  Column 

f  Joint  nuabep 

2  Load  condition  nuaber 

3  Load  in  X-direction 

4  Load  in  Y-direction 

5  Load  in  2-direction 

6  Moaeat  about  X-axis 

7  Moaeat  about  Y-axis 

8  Moaeat  about  Z-axis 

Thi?  series  of  cards  must  be  terminated  by  row  of  alter¬ 
nating  zeroes  and  blanks. 


♦  ♦  ♦ 

ADDS?, HI  or  ADDSF, HI , M2 

This  operation  forms  the  total  stiffness  matrix  naaed  Ml  and 
the  lumped  mass  aatrix  naaed  M2  for  the  structural  system 
from  the  element  stiffness  and  mass  matrices  which  are 
stored  on  low  speed  storage  .  These  matrices  can  be  printed 
with  the  PRINT  operation.  If  M2  is  not  specified,  the  row 
aas s  aatrix  M2  will  net  be  formed. 


ADD  K  ,  HI  ,  M2  ,  M3  ,  N 1 

This  opeyatipn  adds  the  element  stiffness  matrix  named. M2  to 
the  total  stiffness  aatrix  named  Ml,  where  Ml  was  previously 
defined  and  initially  set  to  zero.  M3  is  the  name  of  the 
integer  array  in  which  the  column  nuaber  HI  contains  the  row 
or  column  numbers  in  the  total  stiffness  matrix  where  the 
element  stiffness  terms  are  to  be  added. 


DISPL,M1,M2 


s  operation 


ints  the  displacement 


„oint  Sequence  otder,  where  M2  is  the 
condition  array. 


vector  named  Ml  in 
name  of  the  boundary 
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MEHFBC,M1,M2,M3,H4,H1 

this  operation  aultiplies  the  element  stiffness  aatrix. naaed 
HI  by  the  jo?.nt  displaceaent  aatrix  naaed  M2.  m3  is  the 
na>e  of  the  integer  array  m  which  the  colaan  nuaber  HI 
contains  the  row  nuab^rs  in  the  displaceaent  matrix,  M2, 
which  are  to  be  multiplied  by  the  eleaent  stiffness  (or 
force-displaceaent)  aatrix.  Ml.  The  results  of  this  aulti- 
plication  are  stored  in  the  array  naaed  M4. 


FORCE, Ml, M2, M3  or  FORCE, Ml,  M2 

This  operation  calculates  the  aeaber  forces  for  a  group  of 
eleaent s  m  which 

Ml  is  the  nane  of  the  eleaent  group 
M2  is  the  naae  or  the  displaceaent  aatrix 
M3  is  the  naae  of  the  aatrix  in  which  the  forces  are 
stored  in  the  order  calculated. 

If.  M3  is  not.  specified,  the  .eleaent  forces  will  be,  printed 
only  ana  will  not  Be  retained  in  storage.  For  the  TRUSS 
eleaent  only  the  aeaber  axial  force,  ?,  will  be  calculated 
for  each  aeaber.  For  the  BE&M  eleaent,  eight  forces  will  be 
printed  with  reference  to  the  positive  definition  shown  in 
the  BEAM  operation. 
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3.  si Basis  Aa&iisis  Ssasaiisas 


The  following  operations  were  designed  to  evaluate 
the  dynamic  response  of  structures  subjected  to  arbitrary 
time- dependent  loads.  If  these  operations  are  used  in 
connection  with  the  standard  aatrix  operations  and  the 
structural  analysis  operations,  a  dynaaic  analysis  is  a 
relatively  siaple  procedure.  The  user  has  the  option  of 
using  the  node  superposition  aethod  or  a  direct  step-by-step 
integration  of  the  dynaaic  equations  of  action.  The  user 
may  examine  the  spectra  of  both  input  loading  and  calculated 
displacements.  In  addition,  the  contributions  of  the  indi¬ 
vidual  modes  may  be  evaluated  and  compared. 

The  most  common  and  convenient  form  for  time- 
dependent  data  to  be  specified  is  as  straight  line  segments 
between  given  time  points.  Therefore,  an  operation  which 
generates  values  at  equal  intervals  is  necessary.  Another 
common  characteristic  of  time-varying  loads  on  structures  is 
that  it  is  normally  possible  to  represent  the  loads  at  all 
points  on  the  structure  by  the  product  of  two  matrices,  a 
column  matrix  indicating  the  spatial  distribution  of  loads 
times  a  row  aatrix  which  indicates  the  values  as  a  function 
of  various  times.  If  a  more  complicated  loading  is 
required,  it  is  possible  to  perform  more  analyses,  each 
within  the  restrictions  of  the  program,  then  add  the  results 
of  each  analysis. 

The  following  operations  have  been  added  for  the 
major  purpose  of  performing  dynamic  analysis. 
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POTG, 81, 82,83,81, 82 


This  operation  generates  a 
maes  at  equal  internals 
array  naaed  81.  The  array 
fora: 


?atrix.  naaed  82  which  contains 
the  function  specified  in  the 
81  aust  be  a  2  by  K  array  of  the 


which  nuaerically  represents  a  function  of  the  fora  shown 
below: 


F  (T  ) 


The  tine  internal  t  is  specified  in  the  1  by  1  matrix  naaed 
83.  81  specifies  the  total  number  of  values  to  be  gener¬ 
ated,  and  is  the  nuaber  of  columns  in  82,  If  82  the 
array  82  will  be  a  1  x  8  1  row  matrix  m  which  the  fipst 
value  will  be  f.  If  82  is  not  equal  to  0,  the  array  82  will 
be  a  2  x  81  matrix  of  the  following  fora: 


[82]  - 


ti  t»  ♦  dt  t,  ♦  2dt 
ft  f  (t,  ♦  At)  f  (t,  ♦  2At) 
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•  •  ♦ 

STB  P  ,  Ml  , M2,H3,H4,M5,M6,M7,  M 8,  N1  ,N2 

This  operation  calculates  the  dynamic  response  of  a  struc¬ 
tural  system  using  direct  step-by-step  integration  of  the 
following  linear  aatrix  equation  or  aotion: 

♦  [CKO}  ♦  *  R(t)  =  PF(t) 

She  re:  .  ,  . 

Hi  fs  the  naae  of  the  H  x  N  stiffness  aatrix  K 

H2  is  the  naae  of  the  N  x  H  mass  aatrix  N 

H3  is  the  naae  of  the  N  x  M  damping  aatrix  C 

H4  is  the  naae  of  the  H  x  3  initial  condition  matrix 

a  in  which: 

0  (1,1)  is  a  vector  of  disolacements  0 

0  (1,2)  is  a  vector  of  velocities  0 

D  (1,3)  is  a  vector  of  accelerations  0 

H 5  is  the  naae  of, the  )(  by  N2  matrix  of  calculated 

displacements  in  which  column  i  represents  the 
displacements  at  tiae  i*in*  t 
H6  is  the  naae  of  the  N  x  1  load  distribution  aatrix 
p 

H7  is  the  naae  of  the  1  x  k  row  matrix  representing 
the  load  multipliers  at  equal  time  increments 
?,  where  k  *  Nz/HI 

H8  j.s  the  naae  of  the  1x1  aatrix  containing  t 
HI  is  the  output  interval  fgr  the  displacements 
N2  is  the  total  number  of  displacement  vectors  to 
be  calculated. 


The  total  ti*«  for  which  results  wj.ll  be  calculated  by  this 
operation  is  S1*N2*  t.  This  operation  must  Be  followed  with 
one  data  card  in  free  format  containing  the  following  infor¬ 
mation: 


Item 

2 

3 


Contents 

DELTA 

ALPHA 

THETA 


Different  values  of  delta,  alpha  and  theta  will  allow  the 
user  to  select  different  methods  of  step-by-step  integra¬ 
tion.  The  following  table  lists  some  possibilities: 


DELTA  ALPHA  THETA 


Nevmarks  Average  Acceleration  1/2  1/4  1.0 

Linear  Acceleration  1/2  1/6  1.0 

Wilson's  Theta  Method  (low  damping)  1/2  1/6  1.42 

Wilson's  Theta  Method  (high  damping)  1/2  1/6  2.0 
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EIGEN, 91,92, S3,  N1 


This  operation  solves  the  following  eigenvalue  problem: 


K  « 


9  *  A 


In  which  the  N  xLN,  smtric,  posi*ive-semidefinit§  matrix  K 
Is  pamea  91.  The  matrix  9  is  a  diagonal  matrix  or  nonzero, 
positive  terms  designated  by  93.  The  matrix  93  must  be  a 
row  or  column  matrix  containing  only  the  diagonal  terms  of 
9.  The  eigenvalues,  ,  are  stored  in  matrix  93.  The  eigen¬ 
values  are  ordered  in  numerically  increasing  order  and  the 
eigenvectors,  ,  are  stored  in  the  corresponding  columns  of 
the  matrix  92,  The  number  N1  specifies  the  approximate 
number  of  significant  figures  of  the  eigenvalues.  If  N1  is 
zero  or  blanlc.  4  figure  accuracy  will  be  used, 
accuracy  possible  is  16  figures.  The  use  of 
figure  accuracy  is  not  recommended. 


The  maximum 
more  than  12 


The  program  reduces  the  problem 
by  tne  following  transformation 


to  standard  eigenvalue  form 


where 


K«  a  m  T  K  m 
I  *  a  T  9  m 

in  which 

IB|  *  1  /  /  Hi  i 

The  calculated  mode  shapes,  ,  are  normalized  as  follows: 

=  i  4>tK4>  =  X 

The  program  uses  the  standard  Jacobi  diagonalization  method 
to  solve  for  all  eigenvalues  and  eigenvectors. 
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DYN  AH,M1,M2,M3,  H4,H5,  86,  81 


operation  evaluates  the  .following.  se+  of  uncoupled 
d  order  differential  equations  associated  with  the  Sode 
position  aethod  for  the  dynamic  analysis  of  a  struc- 


This  operati 
second  order 
superposition 
tural  system. 


x,  +  2 A,aj, x,  ♦  ufa,  =J}(t)  i  s  1  to  N  nodes 


H 1  is  the  naae  of  a  row  or  column  matrix  which  contains  the 
N  teras  (frequencies  in  rad/sec)  .  M2  is  the  naae  of  a  row 
or  column  matrix  which  contains  the  8  terms  (ratio  of 
modal  damping  to  critical  damping). 

The  generalized  time-varying  forces  P  (t)  are  not  specified 
directly  But  are  evaluated  from  more  fundamental  informa¬ 
tion.  The  forces  for  all  modes  are  evaluated  at  specific 
times  by  the  program  from  the  following  matrix  equation: 


p  =  p*f 


M3*H4 


In  which  p  is  a  specified  J  :  1  vector  named  S3,  and  f  is  a 
1  x  Ii  row  matrix  which  will  be  generated  from  the  2  By  k 
array  named  B4.  The  array  M4  is  the  same  fora  as  the  input 
array  described  under  the  operation  FUNG.  It  is  not  neces¬ 
sary  to  use  FUNG  before  the  DINAH  operation. 

H5  is  th?  naae  of  the  V  i  11  array  which  contains  the  gener¬ 
alized  displacement  X  (t)  . 

;6  is  the  name  of  the. 1  x  1  array  which  contains  the  time 
cerement  associated  with  the  generalized  displacements. 

N1  is  the  number  of  displacements  to  be  generated. 

The  method  of  integration  used  is  exact  for  straight  line 
segments. 

p£ot,hT7ni 

This  operation  will  prepare  a  printer  plot  of  selective  rows 
of  the  matrix  named  Ml.  N1  is  the  number  of  rows  of  Ml 
which  will  be  plotted  by  this  operation.  This  operation  is 
followed  by  N1  cards  in  f  (1a1,i4)  format  with  the  following 
information: 


Columns  Contain 

1  Plot  symbol 

e  i?  m _  _ i  i _ 


Row  number  to  be  plo 


any  keypunch  symbol 
be  plotted 


The  program  automatically  searches  the  information  to  be 
plotted  for  the  maximum  and  minimum  values.  The  difference 
in  these  numbers  divided  by  120  spaces  is  selected  as  the 
plot  scale. 
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4.  H§£t  itaasfe t  Ofissitjoas 


The  purpose  of  this  series  of  operations  is  to  fora 
the  total  conductivity  and  heat  capacity  matrices  for 
systeas  of  two-  or  three-dimensional  elements,  to  form  the 
flux  vector  and  solve  the  defined  set  of  equations.  Por  two 
dimensional  eleaents  there  are  4  to  9  node  isoparaaetric 
eleaents.  For  three  dimensional  elements  there  are  8  to  21 
node  isoparaatric  elements. 

After  the  creation  of  an  array  containing  the  coor¬ 
dinates  of  the  systea  nodes,  the  specification  of  element 
connectivity,  the  specification  of  material  properties  and 
the  specification  of  constant  temperature  nodes,  the  conduc¬ 
tivity  and  capacity  matrices  and  the  flux  vector  are  formed. 
The  equations  are  solved  by  the  appropriate  equation  solver 
and  the  temperatures  are  printed  in  node  order. 

The  HTXFR  operation  initializes  the  problem.  The 
operation  COOHD  is  used  to  specify  or  generate  nodal  coordi¬ 
nates.  Element  connectivity  is  specified  by  the  EICON  oper¬ 
ation.  Haterial  properties  and  element  boundary  conditions 
are  input  via  the  PROP  operation.  The  CTEHP  operation 
establishes  designated  nodes  as  having  constant  tempera¬ 
tures.  The  equation  profile  for  the  problem  is  generated  by 
the  PROF  operation. 

To  form  the  conductivity  matrix,  the  operation  SYMC 
is  used.  The  heat  capacity  matrix  can  be  approximated  with 
either  a  consistant  (CCAP)  or  lumped  (LCAP)  matrix  formula¬ 
tion.  The  flux  vector  is  formed  with  the  FORM  operation. 

Time  independent  systeas  of  equations  are  solved 
with  the  CALC  operation.  Systems  of  equations  involving  the 
first  derivative  of  temperature  with  respect  to  time  are 
solved  with  the  ODE  series  of  operations.  Nodal  temepra- 
tures  are  printed  in  node  number  order  with  the  PTEKP  opera¬ 
tion. 
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HTXFR 


operation  initializes  the  heat  transfer  problem-  The 
_oj.lowing  information,  entered  m  free  format  must  follow 
this  operation: 


Itea  Contents 

1  Nuaber  of  nodes  (NUMNDP) 

2  Nuaber  of  elements  (NON!) 

3  Nuaber  of  material  sets  (NUMMAT) 

4  Spatial  dimension  (2  or  31  (NDljj 

5  Number  of  unknowns  per  node  (NDF) 

6  Maximum  number  of  nodes  per  element  (NEN) 

The  nuaber  of  unknowns  per  node  will  always  be  1  for  heat 
transfer  problems. 


COORD 

This  operation  creates,  an  array  which  contains  the  coordi¬ 
nates  fot  all  nodes  in  a  heat  transfer  system.  Data  is 
entered  in  free  format  as  follows: 


Itea 

1 

2 

3 

4 

5 

6 


Contents 
Node  nuaber 
X-coordmate 
T-coord inate 
Z -coordinate 
System  type 
Generation  co 


(if  3-D, 
de 


else  oait) 


System  type  refers  to  the  coordinate  system  used  when  input¬ 
ting  data.  All  coordinates  are  converted  to  the  cartesian 
coordinate  system  for  use  by  CAL-NPS. 


System  Type 

V 

3 

4 

5 


System 
Carte sia  n 
Cy  lindncal. 
Cylindrical, 
Cylindrical, 
Spherica  1 


Z-axis  longitudinal 
Y-axis  longitudinal 
X-axis  longitudinal 


The  input  data  is  the  same  as  shown  above  with  the  following 
correspondance: 

Cartesian  Cylindrical  Spherical 

X  ^  p 


If  the  generation  code  is  not. zero  (0),  the. next  card  is  a 

feneration  vector  for  automatic  node  generation.  It  is  m 
ree  format  as  follows: 


Itea  Contents 

1  Node  number  increment 

2  X  increaent 

3  Y  increment 

4  Z  increaent  (if  3-D,  else  omit) 

5  Last  node  number  to  be  generated 

This  operation  aust  be  terminated  by  a  card  with  alternating 
zeroes  ana  blanks. 
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SIC  01 


This  operation .creates  an.  array  which  contains  the  element 
connectivity  data  for  the  elements  in  a  .heat  transfer 
systea.  Data  is  entered  in  free  forest  as  follows: 


Item 


1 


etc 

H-*  1 
H+ 
!♦ 


Contents 
Elesent  nnaber 
lode  1  nnaber 
lode  2  nnaber 
etc 

lode  1  nnaber 
Baterial  set  nnaber 
Generation  code 


If  the  Generation  code  is  not  zero.  (0)  ,  the  next  card  is  a 
generation  vector  of  the  automatic  generation  of  eleaent 
connectivity.  It,  is  entered  in  free  foraat  as  follows: 

Itea  Contents 

1  Eleaent .nnaber  increment 

2  lode  1  lucre ■ ent 

3  lode  2  increment 

etc  etc 

1+1  lode  l  increment 

H+2  Baterial  set  increaent  (nsnally  0) 

H+3  last  eleaent  nnaber  to  be  generated 

This  operation  post  be  terminated  by  a. row  of  alternating 
zeroes  ana  blanks.  The  node  nna Bering  conventions  for 
eleaent  connectivity  are  shown  below. 


2-D  Eleaent 


3-D  Eleaent 
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PROP 


This  operation  inputs  the  material  property  data  for  a  heat 
transfer  system.  The  following  information,  entered  in  free 
format,  must  follow  this  operation: 


Item  Contents 

1  Material  set  number 

2  Element  type  number  (2  for  2-D,  3  for  3-D) 


Additional  information  must  be  provided,  depending  on  the 
type  of  element  being  used.  Note  that  material  property 
data  must  be  in  consistent  units  as  shown  below: 


k 

c 

rho 


BTO/lbm-op 

lbm/ft3 


SI  Units 

_  W/Tff=*C 
kJ/kg-°C 
kg/m3 
R/m3 


2-D  Elements 

This  information  fust  follow  the  PROP  operation  card  an$  its 
two  required  entries.  An  entry  must  be  made  for  each  item. 
If  the  item  is  temperature  dependent,  the  entry  will  be 
ignored. 

Item  Content? 

1  Conduc"  iv^ty  in  the  X-direction 

2  Conductivity  m  the  Y-direction 

3  Specific  Heat 

4  Density 

5  Heat  generation  per  unit  volume 

6  Number  of  integration  points  per  direction 
(1  to  6,  default  is  4) 

7  Geometry  type  {see  below) 

8  Total  number  of  lines  with  specified  boundary 
conditions  in  elements  with  the  same  material 
set  numer  (NLBC)  {see  below) 

9  Temperature  dependence  code  (see  below) 

Geometry  type  is  1  for  plane  geometry  and  2  for  axisymmetry. 

The  temperature  dependence  gode  is  0  if  all  material  proper¬ 
ties  are  constant  and  1  if  any  property  is  temperature 
dependent.  If  the  code  is  1,  the  following  information  (in 
fifee  format)  is  required. 


Item  Content? 

1  Conductivity  j.n  the  X-d3.rect3.on  code 

2  Conductivity  m  the  Y-direction  code 

3  Heat  capacity  (specific  heat'' density)  code 

4  Heat  generated  per  unit  volume  code 

where  0  means  a  constant  property  and  1  means  a  temperature 
dependent  property.  Temperature  dependent  properties  are 
entered  in  the  form  of  a  table.  Tne  tables  are  consecu¬ 
tively  input  for  conductivity  in  the  X-direction,  conduc¬ 
tivity  in  the  Y-direction.  heat  capacity  and  heat  generated 

Ser  unit  volume.  Omit  the  tables  for  which  the  temperature 
ependence  code  is  zero.  Tables  are  input  in  free  format  as 
shown  on  the  following  page. 
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Item 

1 


2^5+1 


Contest  s 

Humber  of  data  pairs  to  be  entered  (This 
should  be  a  single  card) 

Temperature  1 

Property  1  (These  two  entries  should  be  on 
one  card) 

Temperature  2 

Property  2  (These  two  entries  should  be  on 

one  card) 

etc 


Temperature  1! 


one  care 


(These  tvo  entries  should  be  on 


If  any  lines  hatreds  specified  boundary  condition  (HLBC>0),.a 
card  mustbe  submitted  for  each  line.  If  the  same  line  is 
subjected  to  more  than  one  boundary  condition,  a  card  must 
be  used  for  each  one  of  these  conditions.  The  total  number 
of  cards  must  equal  HLBC.  The  information  is  entered  in 
free  format  as  follows: 


free  format  as  follows: 

Item  Contents 

1  Element  number 

2  Boundary  condition  code  (see  below) 

3  Line  code  (see  below) 

4  Property  value  (see  below) 

5  lmblept  temperature 
The  boundary  condition  codes  are: 

1  Plus  ,  , 

2  Convection  (constant  coefficient) 

3  Radiation 

4  Convection  (temperature  dependent  property) 
The  line  codes  are: 

1  S  »  ♦!  line 

-1  S  ■  -1  line 

2  T  »  ♦!  line 

-2  T  *  -1  line 


The  property  values  are: 

Flux  -  Flux  per  unit  area  ,  . 

Convection  -  constant  heat  transfer  coefficient 
(ignored  if  temperature  dependent) 

Badiation  -  Product  of  emissivity  by  Stephan- 
BOltzman  constant 

The  ambient  temperature  is  ignored  for  the  flux  boundary 
condition. 

If  the  boundary  condition  code  is  4  (temperature  dependent 
heat  transfer  coefficient)  .  a  table  must  follow  for  the 
temperature  dependence.  It  is  input  m  free  format  as  is 
shown  on  the  following  page. 
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Contents 

Nuaber  cf  data  pairs  to  be  entered  (This 
should  be  a  single  card) 

Teaperature  1 

Heat  transfer  coeficient  1  (These  two  entries 
should  be  on  one  card) 

Teaperature  2 

Heat  transfer  coeficient  2  (These  two  entries 
should  be  on  one  card) 

0^C 

Teaperature  N  ,  . 

Heat  transfer  coeficient  N  (These  two  entries 
should  be  on  one  card) 


3-D  Eleaents 

This  information  ®ust  follow  the  PROP  operation  card  and  its 
two  required  entries.  An  entry  must  be  aade  for  each  item. 
If  the  item  is  temperature  dependent,  the  entry  will  be 
ignored. 

Item  Contents  .  . . 

1  Conductivity  in  the  X-direct  gon 

2  Conductivity  m  the  T-direction 

3  Conductivity  in  the  Z-direction 

4  Specific  Heat 

5  Density 

6  Heat  generation  per  unit  volume 

7  Humber  of  integration  points  per  direction 
(1  to  6,  default  is  4) 

8  Geometry  type  (see  below) 

9  Total  nuaber  of  surfaces  with  specified 
boundary  conditions  in  elements  with  the  same 
material  set  numer  (HSBC)  (see  below) 

10  Temperature  dependence  code  (see  below) 

Geometry  type  is  1  for  plane  geometry  and  2  for  axisymmetry. 

The  temperature  dependence  ?ode  is  0  if  all  material  proper¬ 
ties  are  constant  and  1  if  any  property  is  temperature 
dependent.  If  the  code  is  1,  the  following  information  (in 
f^ee  format)  is  required. 

Item  Content? 

1  Conductivity  in  the  X-dgrection  code 

2  Conductivity  in  the  v-dgrection  code 

3  Conductivity  m  the  ^-direction  code 

4  Heat  capacity  (specific  heat*dersity)  code 

5  Heat  generated  per  unit  volume  code 

where  0  means  a  constant  property  and  1  means  a  temperature 
dependent  property.  Temperature  dependent  properties  ar.-. 
entered  in  the  fora  of  a  table.  The  tables  are  consecu¬ 
tively  input  for  conductivity  in  the  X-direction.  conduc¬ 
tivity  in  the  X-direction,  conductivity  in  the  Z-airection, 
heat  capacity  and  heat  generated  per  unit  volume.  Omit  the 
tables  for  which  the  temperature  dependence  code  is  zero. 
Tables  are  input  in  free  format  as  is  shown  on  the  following 
page. 
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Item  Contents 

1  lumber  of  data  pairs  to  be  entered  (This 
should  be  a  single  card) 

2  Teeperatore  1 

3  Property  1  (These  two  entries  should  be  on 
one  card) 

4  Temperature 

5  Property  2 
one  card) 

etc  etc 
2*H  Teeperatore 
2*»+l  Property  s 
one  card) 
hare 


(These  two  entries  should  be  on 


(These  two  entries  should  be  on 


The  information 


If  any  surfaces  hare  a  specified  boundary  condition 
(HSBC>0)  *  a  card  must  be  submitted  for  each  surface.  If  the 
same  surface  is  subjected  tc  more  than  one  boundary  condi¬ 
tion.  a  card  must  be  used  for  each  one  of  these  conditions. 
The  total  number  of  cards  must  equal  tfSBC. 
is  entered  in  free  format  a  s  follows: 

Item  Contents 

1  Element  number 

2  Boundary  condition  code  (see  below) 

3  Surface  code  (see  below) 

4  Property  value  (see  below) 

5  Ambient  temperature 
The  boundary  condition  codes  are: 

1  Flux 

2  Convection  (constant  coefficient) 


The 


3  Radiation 

4  Convection 
line  codes  are: 


(temperature  dependent  property) 


1 

-1 

-1 

-I 


R 

R 

S 

S 

T 

T 


The  property^yalues  are; 

Plux  -  Flux  per  unit  area 

Convection  -  constant , heat  transfer  coefficient 
(ignored  if  temperature  dependent) 

Radiation  -  Product  of  emissivity  by  Stephan- 
Boltzman  constant 

The  ambient  temperature  is  ignored  for  the  flux  boundary 
condition. 

If  the  boundary  condition  code  is  4  (temperature  dependent 
heat  transfer  coefficient)  .  a  table  must  follow  for  the 
temperature  dependence.  It  is  input  in  free  format  as  is 
shown  on  the  following  page. 
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Item  Contents 

1  Number  of  data  pairs  to  be  entered  (This 
should  be  a  single  card) 

2  Temperature  1 

3  Hear  transfer  coeficient  1  (These  two  entries 

should  be  on  one  card) 

4  Temperature  2 

5  Heat  transfer  coeficient  2  (These  two  entries 

should  be  on  one  card) 

etc  etc 

2*N  Temperature  N 

2*N+1  Heat  transfer  coeficient  N  (These  two  entries 

should  be  on  one  card) 


CTEHP 

This  operation  inputs  .  constant  temperature  boundary 
restraint  data.  If  you  have  nodes  with  constant  tempera¬ 
ture,  you  must  enter  the  temperatures  for  those  nodes.  For 
nodes  with  no  temperature  restrictions,  no  entries  should  be 
made.  Automatic  generation  capability  is  built  into  the 
operation.  Data  is  entered  in  free  format  as  follows: 

Item  Contents 

1  Initial  node 

2  Last  node 

3  Node  increment 

4  Temperature 

These  entries  may  be  repeat ed. until  all  the  constant  temper¬ 
ature  nodes  are  entered.  This  operation  must  be  terminated 
by  a  row  of  alternating  zeroes  and  blanks. 

This  operation  establishes  the  profile  of  the  equations  for 
solution  of  the  heat  transfer  .problem.  After  the  issuance 
of  this  command,  the  problem  is  set  and  you  may  not  change 
the  node  numbers  with  restrained  boundary  conditions 
(constant  temperatures).  The  values  of  the  restrained  temp¬ 
eratures  may  be  changed. 


SYHC 

This  operation  forms  the  symmetric  conductance  matrix  for 
heat  transfer  problems. 


OSTNC 

This  operation  forms  an  unsymmetric  conductance  matrix. 

rcxp - 

This  operation  forms  a  lumped  capacitance  approximation 
matrix  for  heat  transfer  problems. 
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CCAP 

This  operation  forms  a  consistent  capacitance  approximation 
aatrix  for  Heat  transfer  problems. 


FOBS 


CALC 

I  his  operation  solves  time  independent  heat  transfer  prob- 
ens  for  temperature  and  updates  the  temperature  matrix. 


0DB,M1 

This  operation. solves  the.  first  order  ordinary  differential 
equations  arising  from  time  dependent  heat  transfer  prob¬ 
ieas.  If  Ml  is  '*IBJT,,,the  initial  conditions  are  estab¬ 
lished.  This  operation  must  be  followed  by  the  following 
data  in  free  format: 

Item  Contents 

1  Integration  paraater  theta  for  two  point 
scheme  (default  =  2/3) 

2  Integration  parameter  gamma  for  three  point 
scheme  (default  =  1.5) 

3  Integration  parameter  beta  for  three  point 
scheme  (default  *  0.8) 

4  Maximum  temperature  difference  for  time  step 
adjustment 

5  Minimum  temperature  difference  for  time  step 
ad justm ent 

Default  values  are  obtained  by  entering  zeros  for  items  1r  2 
and  3. 


If 

If 


Ml 

Ml 


|>s  "LIME", 
[s  "QUAD"; 


the  two  point  integration  scheme  is  used, 
the  three  point  integration  scheme  is  used. 


Some  suggested  values  for  theta  for  the  two-point  scheme  as 
given  in  Reference  4  are: 


Ctank-Nicolson 

Zienkiewicz 

Bettencourt 

Liniger 


THETA 

2/3 

^78 


Some  suggested  values  for  beta  and  gamma  for  the  three  point 
scheme  as  given  in  Reference  4  are: 


Lees 

ssir, 

Zienkiewicz 

Bettencourt 


Beta 

1/3 

3/4 

.646 

4/5 

9/10 


Gamma 

1/2 

1.0 

1.2184 

3/2 

3/2 
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PTE  HP 


This  operation  prints  the  nodal  tea pera tares  of  a  hear 
transfer  system  in  node  nuaber  order. 


TOL ,  H 1 

This  operation  sets  the  solution  conyergence  tolerance  to 
the  value  round  m  the  1x1  matrix  named  Hi. 


cony 

This  operation  performs  a  temperature  convergence  test  on  a 
heat  transfer  system.  If  this  operation  is  used  inside  a 
loop  (LOOP  operation)  and  the  test  shows  convergence, 
looping  will  be  terainated. 


DTI  H,  HI 

This  operation  sets  the  tine  increment  fo£  integration  in  a 
heat  transfer  system  to  the  value  found  m  the  1x1  matrix 

Hie 


ADTIH 

This  operation  advances  the. time  in  a  heat  transfer  problem 
by  one  tine  step.  The  time  step  Is  input  with  the  DTIH 
operation. 


EIGV 

This  operation  conputes  the  dominant  eigenvalue  and  eigen¬ 
vector  of  the  current  heat  transfer  conductance  matrix. 


PRO  HPT 

This  operation  permits  the  prompts  for  user  input  to  be 
suppressed  without  loss  or  other  output.  A  second  PROHPT 
will  restore  the  prompts  unless  the  print  output  is 
suppressed  by  the  NO  coamand. 


5.  graphic!  Operations 


CAL-HPS  has  a  Halted  graphics,  capability  for  users 
with  PLOT- 10  coapa tibia  terminals  or  IBB  3277  Dual  Screen 
terminals.  Two-  or  three- diaensional  aeshes  nay  be  viewed 
in  the  X-T,  T-Z  or  X-Z  planes. 

The  graphics  operations  are  initialized  with  the 
GBAPH  operation  by  which  the  user  specifies  which  graphics 
capable  terminal  is  being  used.  The  PLHX  operation  plots 
heat  transfer  aeshes  and  the  PLST  operation  plots  structural 
aeshes.  The  user  aay  title  the  plots  using  the  TITLE  opera¬ 
tion. 


i  initializes  the  graphics . package.  ,  Presently 
terminals  users  aay  utilize  tor  the  graphics 
X  yes  must  be  input  immediately  following  th^s 


GBAPB 

This  operation  initializes  the  graphics 
there  £re  two  " -  -  — 

Operations.  a  jna  ■  wk  U3  l  x.ui 3  uxa  ;  i _ _ ^  _ 

card  to  signify  that  one  of  the  two  types  of  terminals  is 
being  used.  Then  the  terminal  code  must  be  entered  in  free 
format:  ...  J 

Terminal  Code  Teramal 

1  PLOT-10  Compatible  Terminal 

2  IBB  327  7  Dual  Screen 


PLHX, SI 


This  operation  plots  2-D  and  3-D  heat  transfer  analysis 
aeshes.  If  the  mesh  is  3-D,  the  user  aay  specify  the  plane 
cf  view: 

SI  *  1  :  X-I  plane 

SI  »  2  :  I-Z  plane 

Si  *  3  :  X-Z  plane 

The  default  is  the  I-T  plane, 

2-D  aeshes. 


hence  SI  is  not  required  for 


PLST,H1,H2#H1 


This  operation  plots  2-0  aad  3-D  structural  analysis  meshes. 
Hi  is  the  nods  coordinate  aatriz  created  by  the  BODES  opera¬ 
tion.  82  is  the  ala aent  connectivity  aatriz.  loaded  with 
the  LOIDI  operation.  The  row  dimension  is  the  number  of 
elaaents  and  the  col  nan  diaension  is  the  aaziaua  nnaber  of 
nodes  per  elesent  plus  one.  82  contains  tha  following 
information: 

H2(K,  1)  *  Humber  of  nodes  in  element  n amber  K 


H2(K,  1) 
etc. 


■  number 
*  Bode  1 

■  Bode  2 

etc. 


82 (K, B*1)  *  Bode  H 
HI  specifies  the  plane  of  view: 
HI  *  1  :  X-Y  plane 

HI  *  2  :  T-Z  plane 

HI  »  3  :  X-Z  plane 

If  tha  structure  contains  two 
(PL ABE  operation),  the  connect 
tion  shown  below: 


e  contains  two  dimensional  aeabrane  elaaents 
a) ,  the  connectivity  east  follow  the  conven- 


N.  H  Ji 


TITLE  ,81 


s  operation  allows  the 
teen  character  lines  to 


the  PLHX  or  PLST  op 
upper  right  hand  a 
area.  .  This  operat 


!ws  the  user  to  input  HI  (up  to  three) 
nes  to  label  plots  generated  with  either 
erations.  The  label  will  appear  in  the 
rea  of  the  screen  outside  the  plotting 


operation.  1 
command.  The 


operation  aust  immediately  preceed  the  plotting 
The  title  may  be  changed  by  reissuing  the  TITLE 
default  value  of  HI  is  one. 
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6 .  looping  operations 

CAL- UPS  has  a  five  level  looping  ability.  The  first 
operation  is  LOOP  and  the  last  operation  is  NEXT. 
Operations  within  CAL-NPS  are  noraally  executed  as  they  are 
encountered.  If  the  operation  requires  data,  the  data  cards 
follow  the  operation  card.  In  the  case  of  looping,  however, 
all  operation  cards  are  stored  within  the  computer  before 
they  are  executed.  If  operations  within  the  loop  require 
data,  the  data  cards  aust  be  supplied  in  the  order  required 
after  the  last  NEXT  operation.  If  an  error  is  encountered 
while  executing  in  a  loop,  the  entire  matrix  of  loop 
cosaauds  is  deleted  and  the  user  is  given  the  opportunity  to 
try  again.  Hatrices  that  have  been  aodified  by  operations 
successfully  coapleted  while  in  the  loop  remain  aodified. 
After  all  loops  are  executed  the  computer  storage  required 
for  these  operations  is  automatically  released  by  the 
program. 
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LOOP, HI 

HI  is  the  number  of  tiaes  the  loop  is  to  be  executed. 
Associated  with  each  LOOP  operation  there  must  be  a  corre¬ 
sponding  NEXT  operation  which  signifies  the  end  of  the  loop 
and  the  return  of  the  control  to  the  beginning  of  the  loop. 
The  following  is  a  possible  series  of  looping  operations. 


LOOP,  5-< 


LOOP,  2* 


NEXT 


LOOP.  4 


HEXT 

NEXT 


vl 


Second  level  loop 
executed  twice 
(total  cf  10  times) 


Second  level  loop 
executed  4  times 
(total  20  times) 


Pirst  level  loop 
executed  5  tiaes 


HEXT, HI  or  NEXT 

The  operation  NEXT  signifies  the  end  of  a  loop,  .  It,  is 
apparent  which  LOOP  ana  NEXT  cards  are  associated  if  there 
are  an  egual  number  of  each.  The  operation  NEXT, Ml  will 
cause  the  loop  to  terminate  if  the  first  term  in  the  matrix 
HI  is  negative. 


SKIP,  HI,  N1 


This  operation  will  cause  the. skip  of 
if  the  first  term  in  the  matrix  naned 


the. next  N 1 . operations 
Ml  is  negative. 


7.  Osqr  Defined  Operat ions 


USER  A  and  USERS 

These  names  are  reserved  for  operations  to  be  defined  and 
programmed  by  the  user.  In  order  to  program  these  opera¬ 
tions  it  is  necessary  to  understand  the  internal  organiza¬ 
tion  of  CAL-NPS.  Chapter  III  of  Reference  2  contains  the 
details. 


B.  LIEGE  PROBLEMS 


CA1-NPS  is  designed  as  an  edacational  tool.  It  does  not 
take  advantage  of  banding  and  symmetry  in  matrix  storage, 
except  in  the  heat  transfer  operations.  Larger  probleas  can 
be  solved  by  increasing  the  diaension  of  the  L  array,  but  a 
general  purpose  program  that  aakes  aaxiaua  advantage  of 
out-of-core  storage  and  takes  advantage  of  banding  and 
syaaetry  for  in-core  matrix  storage  is  probably  a  better 
choice.  With  the  above  disclaiaer,  to  increase  problem  size 
capability,  increase  the  diaension  of  the  L  array  and  change 
the  value  of  MAX  to* the  nev  diaension  size  in  the  following: 

C - MAIN  PROGS  AH 

C - SET  PROGRAM  CAPACITY 

COMMON  HTOT,NDP,  L (100000) 

MTOT  »  100000 
NDP  »  2 
CALL  SBTIME 
CALL  CALI 
STOP 
END 

With  the  dimension  of  the  L  array  as  above,  the  program 
currently  executes  in  1024K  bytes  for  CP/CHS.  The  region 
necessary  for  execution  will  increase  about  eight  tiaes  the 
increase  in  the  L  array. 
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APPENDIX  B 


SAMPLE  DATA  FILE 

This  is  a  saaple  data  file  (FILE  FT04F001)  for  the 
hollow  cylinder  with  circnaferential  heating  strips  problea. 


1 


2  1  4 
COOED 

I  . 1666667  5. 

II  -.0013019 
1666667  4. 
-.0013019 
1666667  3. 
-.0013019 
1666667  2. 
-.0013019 
1666667  2. 
-.0013019 
1666667  2. 
-.0013019 
1666667  2. 
-.0013019 
1666667  1. 
-.0013019 
1666667  1. 
-.0013019 
.1666667  0 
-.0013019 
.1666667  0 
-.0013019 
.161458  5. 
-.003906  0 
.161458  4. 
-.003906  0 
.161458  3. 
-.003906  0 
.161458  2. 
-.  003906  0 
.161458  2. 
-.  003906  0 
.161458  2. 


2 

11 

3 
11 

4 
11 

5 

11 

6 
11 

7 
11 

8 
11 

9 
11 

10 
11 
11 
11 

34 
11 

35 
11 

36 
11 

37 
11 

38 
11 

39 
11 

40 
11 

41 
11 

42 
11 

43 
11 

44 
11 


*.  003906  0 
.161458  2. 
-.003  90  6  0 
.161458  1. 
-.003906  0 
.161458  1. 
-.003906  0 
.161458  0. 
-.003906  0 
.161458  0 
-.003906  0 


819103 
0  23 
252421 
0  24 
46  9080 
0  25 
909551 
0  26 
461928 
0  27 
238116 
0  28 
014305 

?9  0493 
0  30 
342870 
0  31 
.78334 
0  32 
2  1 
0  33 
819103 

56 

252421 

57 

469080 

58 

909551 

59 

461928 

2!8l16 

61 

014305 

62 

79  0493 
63 

342870 

713  341 

65 
2  1 

66 


2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
1  2  1 

2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
2  1 
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67  .1*7135  5.819103 
11  -.0110675  0  89 

68  .1*7135  4.252421 
11  -.0110675  0  90 

69  .147135  3.469080 
11  -.0110675  0  91 

70  .  147135  2.909551 

71  7i471§i72.46ll28 
11  -.0110675  0  93 

72  .147  135  2.23  8116 
11  -.01106  7  5  0  94 

73  .  147135  2.  014305 

74  ri471?|71.?90493 
11  -.0110675  0  96 

75  .147  135  1.34  2  870 
11  -.0110675  0  97 

76  .147135  0.783341 
11  -.0110675  0  98 

77  .147135  0  2  1 
11  -.0110675  0  99 
000000000 
EICON 

I  12  13  2  1  1  1 
1  1  1  1  1  0  10 

II  23  24  13  12  1  1 
11  11  1  0  20 

21  34  35  24  23  1  1 
1  1  1  1  1  0  30 
31  45  46  35  34  1  1 
1  1  1  1  1  0  40 
41  56  57  46  45  1  1 

I  1  1  1  1  0  50 

51  67  68  57  56  1  1 
1  1  1  1  1  0  60 
61  78  79  68  67  1  1 

II  1  1  1  0  70 

71  89  90  79  78  1  1 
1  1  1  1  1  0  80 
0000000000 
CTBHP 
1  6  1  160 
89  99  1  60 
0  0  0  0  0  0 
PROP 
1  2 

.1  .1  .6  70 


0  0 


0  0  0 


4  1 

H 
?? 


0 

10  60 
10  60 
10  60 
10  60 


10  2  2  10  60 

PROP 

SIBC 

FOHH 

CALC 

PTBHP 

READ,  5 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


0  0  0  0  0  0 


0  0  0  0 
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APPENDIX  C 


r 


EXAMPLE  TERMINAL  SESSION 

The  following  terminal  session  was  recorded  using  the 
data  file  in  Appendix  B. 


ENTER  TERMINAL  CODE: 

3  *  Any  Alpha  Numeric  Terminal  (NO  GRAPHICS) 

3 

IBM  3278  TERMINAL 
CP  TERMINAL  LINESIZE  132 

GLOBAL  TXTLIB  CMSLIE  FORTMO D2  MOD2EEH  NONIMSL  IMSLSP 
PILED EF  01  DISK  FILE  01  (RECFM  VS  LRECL  3408  BLKSIZE  3412) 
FILEDEF  08  DISK  FILE  FT08F0  01  (RECFM  FBA  LRECL  132  BLKSIZE 
FILEDEF  04  DISK 

FILEDEF  13  DISK  CAL  TEST  (RECFM  FBA  LRECL  132  BLKSIZE  132) 
FILEDEF  50  DISK  HLP  CAL  C 

LOAD  CALO  CAL  FRTCMX  GROOP1  GR00P2  GROUP3  GROOP4  GR  (CLEAR 
EXECOTION  BEGINS... 

-  0. 0  SECONDS 

START 

**START 


-  0.003  SECONDS 

PROMPT 

>w»prohpt 

-  0.00  3  SECONDS 

READ.  4 
#*REiDf  4 

- 1 -  0.003  SECONDS 


*#HTXFR 


99 

ROMS 

1 

COLOM  NS 

1 

ROMS 

8 

COL DM  NS 

2 

HOMS 

4 

COLOMNS 

4 

ROMS 

1 

COLOM NS 

1 

ROMS 

4 

COLOMNS 

4 

ROMS 

1 

COLOMNS 

4 

ROMS 

4 

COLOMNS 

7 

ROMS 

1 

COLOM  NS 

10 

ROMS 

1 

COLOMNS 

1 

ROMS 

99 

COLOMNS 

2 

ROMS 

99 

COLOM  NS 

5 

ROMS 

80 

COLOM  NS 

1 

ROMS 

99 

COLOMNS 

99 

ROMS 

1 

COLOMNS 

■;9 

ROMS 

1 

COLOMNS 

NUMBER  OP  NODAL  POINTS 
NUMBER  OF  ELEMENTS  « 

NUMBER  OF  MATERIAL  SETS  * 

DIMENSION  OF  COORDINATE  SPACE= 
DEGREES  OF  FREEDOM/NODB 
NODES  PER  ELEMENT  (MAXIMUM)  = 


••COORD 


99 

80 

1 

2 

1 

4 

003  SECONDS 


X 

l 

0.166 

0.0  17 

2 

0.  166 

0.012 

3 

0.  166 

0.010 

4 

0.  166 

0.008 

5 

0.  167 

0.007 

6 

0.  167 

0.007 

7 

0.  167 

0.006 

8 

0.  167 

0.005 

9 

0.  167 

0.004 

10 

0.  167 

0.002 

11 

0.  167 

0.0 

12 

0.165 

0.017 

13 

0.  165 

0.012 

14 

0.  165 

0.010 

15 

0.  165 

0.008 

16 

0.  165 

0.007 

17 

0.  165 

0.006 

18 

0.  165 

0.006 

19 

0.  165 

0.005 

20 

0.  165 

0.004 

21 

0.165 

0.002 

22 

0.165 

0.0 

23 

0.  163 

0.017 

24 

0.164 

0.012 

25 

0.  164 

0.010 

26 

0.  164 

0.008 

27 

0.  164 

0.007 

28 

0.  164 

0.006 

29 

0.164 

0.006 

30 

0.164 

0.005 

31 

0.  164 

0.004 

32 

0.164 

0.002 

33 

0.  164 

0.0 

34 

0.  161 

0.016 

35 

0.  161 

0.012 

36 

0.  161 

0.010 

37 

0.  161 

0.008 

38 

0.  16  1 

0.007 

39 

0.  161 

0.006 

40 

0.  161 

0.006 

41 

0.  16  1 

0.005 

42 

0.  16  1 

0.004 

43 

0.  161 

0.002 

44 

0.  161 

0.0 

45 

0.157 

0.016 

46 

0.  157 

0.012 

47 

0.  157 

0.010 

48 

0.  157 

0.008 

49 

0.157 

0.007 

50 

0.  157 

0.006 

51 

0.  157 

0.006 

H 

0.  157 

0.005 

0.  158 

0.004 

54 

0.  158 

0.002 

55 

0.  158 

0.0 

56 

0.  153 

0.016 

57 

0. 153 

0.011 

58 

0.  153 

0.009 

81 


59 

0.  153 

60 

0.  154 

61 

0.  154 

62 

0.  154 

63 

0.  154 

64 

0.  154 

65 

0.  154 

66 

0.  154 

67 

0.  146 

68 

0.  147 

69 

0.  147 

70 

0.  147 

71 

0.  147 

72 

0.  147 

73 

Q.  147 

74 

0.  147 

75 

0.  147 

76 

0.  147 

77 

0.  147 

78 

0.  135 

79 

0.  136 

80 

0.  136 

81 

0.136 

82 

0.136 

83 

0.  136 

84 

0.  136 

85 

0.  136 

36 

0.  136 

87 

0.  136 

88 

0.  136 

89 

0.  124 

90 

0.  125 

91 

0.125 

92 

0.  125 

93 

0.125 

94 

0.  125 

95 

0.  12  5 

96 

0.  125 

97 

0.  125 

98 

0.  125 

99 

0.  125 

**ELCON 

1 

12 

13 

2 

1 

1 

2 

13 

14 

3 

2 

1 

3 

14 

15 

4 

3 

1 

4 

15 

16 

5 

4 

1 

5 

16 

17 

6 

5 

1 

6 

17 

18 

7 

6 

1 

7 

18 

19 

8 

7 

1 

8 

19 

20 

9 

8 

1 

9 

20 

21 

10 

9 

1 

10 

21 

22 

11 

10 

1 

1  1 

23 

24 

13 

12 

1 

12 

24 

25 

14 

13 

1 

13 

25 

26 

15 

14 

1 

14 

26 

27 

16 

15 

1 

15 

27 

28 

17 

16 

1 

16 

28 

29 

18 

17 

1 

17 

29 

30 

19 

18 

1 

18 

30 

31 

20 

19 

1 

19 

31 

32 

21 

20 

1 

20 

32 

33 

22 

21 

1 

21 

34 

35 

24 

23 

1 

22 

35 

36 

25 

24 

1 

23 

36 

37 

26 

25 

1 

0.008 

0.007 


0.005 

0.004 

0.002 

0.0 

0.015 

0.011 

0.009 

0.007 

0.006 

0.006 

8:881 
0.003 
0.002 
0.0 
0.014 
0.010 
0.008 
0.007 
'  .006 
,005 
,005 
,004 
.003 
.002 
.0 
0.013 
0.009 
0.008 
0.006 
0.005 
0.005 
0.004 
0.004 
0.003 
0.002 
0.0 


0. 

0. 

0. 

0. 

0. 

0. 

0. 


0.  083  SECONDS 
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24 

37 

38 

27 

26 

1 

25 

38 

39 

28 

27 

1 

26 

39 

40 

29 

28 

1 

27 

40 

41 

30 

29 

1 

28 

41 

42 

31 

30 

1 

29 

42 

43 

32 

31 

1 

30 

43 

44 

33 

32 

1 

3  1 

45 

46 

35 

34 

1 

32 

46 

47 

36 

35 

1 

33 

47 

48 

37 

36 

1 

1 4 

48 

49 

38 

37 

1 

35 

49 

50 

39 

38 

1 

36 

50 

51 

40 

39 

1 

37 

51 

52 

41 

40 

1 

38 

52 

53 

42 

41 

1 

39 

53 

54 

43 

42 

1 

40 

54 

55 

44 

43 

1 

4  1 

56 

57 

46 

45 

1 

42 

57 

58 

47 

46 

1 

43 

58 

59 

48 

47 

1 

44 

59 

60 

49 

48 

1 

45 

60 

61 

50 

49 

1 

46 

61 

62 

51 

50 

1 

47 

62 

63 

52 

51 

1 

48 

63 

64 

53 

52 

1 

49 

64 

65 

54 

53 

1 

50 

65 

66 

55 

54 

1 

51 

67 

68 

57 

56 

1 

52 

68 

69 

58 

57 

1 

53 

69 

70 

59 

58 

1 

54 

70 

71 

60 

59 

1 

55 

71 

72 

61 

60 

1 

56 

72 

73 

62 

61 

1 

57 

73 

74 

63 

62 

1 

58 

74 

75 

64 

63 

1 

59 

75 

76 

65 

64 

1 

60 

76 

77 

66 

65 

1 

6  1 

78 

79 

68 

67 

1 

62 

79 

80 

69 

68 

1 

63 

80 

81 

70 

69 

1 

64 

81 

82 

71 

70 

1 

65 

82 

83 

72 

71 

1 

66 

83 

84 

73 

72 

1 

67 

84 

85 

74 

73 

1 

68 

85 

86 

75 

74 

1 

69 

86 

87 

76 

75 

1 

70 

87 

88 

77 

76 

1 

71 

89 

90 

79 

78 

1 

72 

90 

91 

80 

79 

1 

73 

91 

92 

81 

80 

1 

74 

92 

93 

82 

81 

1 

75 

93 

94 

83 

82 

1 

76 

94 

95 

84 

83 

1 

77 

95 

96 

85 

84 

1 

78 

96 

97 

86 

85 

1 

79 

97 

98 

87 

86 

1 

80 

98 

99 

88 

87 

1 

**CTEHP 

NODE 

1 

2 

3 

4 

5 

6 


TEBPEBgTggE 

160l00 

160.00 

160.00 

160.00 

160.00 


0.037  SECONDS 
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NODE 

TEMPERATURE 

89 

6  0.000 

90 

60.000 

91 

60.000 

92 

6  0.000 

93 

6  0.000 

94 

60.000 

95 

60.000 

96 

60.000 

97 

60.000 

98 

60.000 

99 

60.000 

«*PROP 

0.007  SECONDS 


MATERIAL  SET  1  FOR  ELEMENT  TYPE  2 
DEGREE  OF  FREEDOM  ASSIGNMENTS 


HEAT  CONDUCTION  ELEMENT 


KX  = 


LOCAL 

NUMBER 

1 


GLOBAL 

NUMBER 

1 


.  1 

.  6001 
70. 00000 
.0 


00 


KY  * 


5 

3 

5 


COND DCTIVITY 
SPECIFIC  HEAT 
DENSITY 

HEAT  GENER/UNIT  VOL 

4  GAUSS  PT S/DIR 

5  LINES  BITH  SPECIFIED  BOUNDARY  CONDITIONS 
PLANE  ANALYSIS 


.  1000 


ROBS 

ROBS 

ROBS 


3  COLUMNS 
9  COLUMNS 
2  COLUMNS 


LINE  B.C. 

ELEg  B.5.  LIJE  PR^E^JALUE  Jglgg  |JJ"« 

7  2  2  10.  00000  60.00000 

8  2  2  10.  00000  60.00000 

9  2  2  10.  00000  60.00000 

10  2  2  10.  00000  60.00000 

-  0.027  SECONDS 

’M'PROF 

8  2  ROBS  1  COLUMNS 

-  0.017  SECONDS 

•fcSYMC 

94  5  ROBS  1  COLUMNS 

-  0.406  SECONDS 

»*FORM 

FLUX  CONVERGENCE  TEST 

RNHAX  »  60.927  RN  =  60.927 

-  0.716  SECONDS 

*«*CALC 

ENERGY  (DB*A<*DR)  =  0.  2 2491 7705 2D  +  0 5 

- -I -  0.030  SECONDS 

*<«PTEHP 

NODAL  TEMPERATURES  TIME  0.0 


NODE 

1 

2 

3 

4 

5 

6 
7 


0.  i255Sd+03 

0.  16  00  0D+03 
0. 16000D*Q3 
0.  160000*03 
0.  16  00 OD* 03 
0.  160000*03 
0.  139700*03 
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8 

0.  12995D+  03 

9 

0.  119740*03 

10 

0.  11  36  2D*  03 

11 

0. 11094D+03 

12 

0.  15581D+03 

13 

0.  1552  6D*03 

1« 

Q.  1541  9D+03 

15 

0.1523 3D*03 

16 

0.  14821D+03 

17 

0. 143170*03 

18 

0.  139920*03 

19 

0.  13  43  4D+Q3 

20 

0.  1256  8D+03 

21 

0.  119570*03 

22 

0.  116800*03 

23 

0.1516  30*03 

24 

0.  150640*03 

25 

0. 148800*03 

26 

0.  14  597D*  03 

27 

0.14 199D+03 

28 

0.  13  96  2D*  03 

29 

0.  1366  4D*  03 

30 

0.  13384D*  03 

31 

0.  12838D+03 

32 

0.  1234  90*03 

33 

0.12  10  4D+03 

34 

0.14  35  5D+03 

35 

0.142140*03 

36 

0.  13999D*  03 

37 

0. 13  76  3D*03 

38 

0.  1353  0D*  03 

39 

0. 13400D+03 

40 

0.  13274D+03 

41 

0.  13142D+03 

42 

0.  128970*03 

43 

0. 12656D+03 

44 

0. 1252  0D+03 

45 

0.132570*03 

46 

0.  131480*03 

47 

0.  13016D*03 

48 

0.  128980*03 

49 

0.  12795D+03 

50 

0.  1274  4D+03 

NODAL  TEMPERATURES  TIME 


NODE 

51 

TEMP 

0. 12693D*03 

52 

0.  1264  4D+03 

53 

0.12556D+03 

54 

0.124730* 03 

55 

0. 124270*03 

56 

0. 1230 OD* 03 

57 

0.  1224  4D*03 

58 

0. 1218  ID*  03 

59 

0. 12 127D*03 

60 

0.  12084D*  03 

61 

0.  1206  20*03 

62 

0.  1204  ID* 03 

63 

0.  1202  ID* 03 

64 

0. 1198 7D* 03 

65 

0. 11955D+03 

66 

Q.  11  93  7D+03 

67 

0.  10875D+03 

68 

0.  1086  20*03 

69 

0.  1084  8D+03 

70 

0.  10  83  6D+03 

!•  M 
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